< 



43 



en 



LU-ITP 2010/003 

IPhT-T10/92 

LAPTH-026/10 



From correlation functions to scattering amplitudes 

Burkhard Eden*, Gregory P. Korchemsky^, Emery Sokatchev* 



* Institut fur theoretische Physik, Universitdt Leipzig 
Postfach 100920, D- 04009 Leipzig, Germany 

t Institut de Physique Theorique 1 , CEA Saclay, 
91191 Gif-sur-Yvette Cedex, France 

* LAPTH 2 , Unwersite de Savoie, CNRS, 
B.P. 110, F-74941 Annecy-le-Vieux, France 



> 

CA ■ Abstract 



We study the correlators of half-BPS protected operators in M = 4 super- Yang-Mills theory, in 
the limit where the positions of the adjacent operators become light-like separated. We compute 
the loop corrections by means of Lagrangian insertions. The divergences resulting from the 
light-cone limit are regularized by changing the dimension of the integration measure over the 
insertion points. Switching from coordinates to dual momenta, we show that the logarithm of the 
correlator is identical with twice the logarithm of the matching MHV gluon scattering amplitude. 
c$ ' We present a number of examples of this new relation, at one and two loops. 
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1 Introduction 

One of the most remarkable manifestations of the AdS/CFT correspondence [1] in the recent 
years was the duality between planar gluon scattering amplitudes and light-like polygonal Wilson 
loops. It was first proposed at strong coupling [2] and soon afterwards also observed at weak 
coupling, first at one loop [3, 4], followed by extensive two-loop tests [5, 6, 7, 8]. This duality 
can be formulated as follows: 

In (AJA^) = In (W[C n }) + 0(1/N C ) + 0(e) . (1.1) 

Here A n is the all-order n— gluon MHV scattering amplitude depending on the particle light-like 
four- momenta p,i (with p\ = and X^=i Pi = 0)> an< ^ ^ n ls ^ ne tree- level amplitude. An essential 
step in establishing the relation (1.1) is the so-called T-duality transformation from momenta to 
dual coordinates: 

Pi %i *Ej+l = *Ei,i+l i -£j j+l U , 2Jj+n = %i ■ \1 - ^/ 

The Wilson loop W[C„] is defined on a closed polygonal contour C n in the dual space, with 
cusps at points Xi and with light-like sides [xj,x i+ i]. The duality (1.1) holds in planar J\f = 4 
super- Yang-Mills theory (SYM) and up to terms vanishing if the suitably identified infrared (for 
amplitudes) and ultraviolet (for Wilson loops) regulators e — > 0. 1 

An important ingredient in this duality is the notion of dual conformal symmetry. This 
is the natural symmetry of the light-like Wilson loop, becoming anomalous due to the cusp 
singularities [10, 11, 12]. By virtue of the duality (1.1) it is extended to a symmetry of the planar 
scattering amplitudes of dynamical origin. The first evidence for this new symmetry came from 
the study [13] of the loop momentum integrals appearing in the four-gluon amplitudes up to four 
(or even five) loops [14, 15, 16, 17]. Once rewritten in dual space according to (1.2), they become 
pseudo-conformal (the infrared regulator breaks the symmetry). This dual conformal symmetry, 
or rather its anomalous version [6, 18], was instrumental in explaining the so-called BDS ansatz 
for MHV amplitudes [15]. 

In the present paper we provide evidence for another duality relation in the M = 4 SYM 
theory, this time between the MHV gluon amplitudes and the correlation functions of gauge 
invariant composite operators on the light cone. The operators we consider belong to the class 
of half-BPS (or "short") scalar operators. They are of the type O^ = Tr(0 fc ), made of the 
six real scalars <fi of the M = 4 SYM theory. They carry R-symmetry SU(4) Dynkin labels 
[0, A;, 0] and transform as chiral primaries under the super conformal symmetry PSU(2,2\4) of 
the M = 4 theory, with fixed conformal dimension d = k. In perturbation theory, such operators 
do not undergo renormalization and are thus protected to all orders. The best known example is 
the simplest, bilinear (k = 2) operator, belonging to the so-called stress-tensor super conformal 
multiplet. In the context of the AdS/CFT correspondence, these operators are dual to massive 
Kaluza-Klein modes in the compactification of type IIB supergravity on an AdSs x S 5 background. 

The correlators of half-BPS operators have been the subject of numerous studies. Not only the 
conformal dimension of the operators, but also their two- and three-point correlation functions 
are protected [19]. The first non-trivial quantum corrections appear in the four-point correlators 



1 Quite remarkably, the duality between planar amplitudes and light-like Wilson loops also holds in gauge 
theories with less or no supersymmctry, including QCD. However, in distinction with J\f = 4 SYM, there the 
relation (1.1) is satisfied in the high-energy (Regge) limit only [9, 3]. 



of protected bilinear operators, which have been computed up to two loops in [20, 21, 22, 23, 24]. 
The knowledge of these quantum corrections allowed one to extract the spectrum of anomalous 
dimensions of the Konishi operator [25] , and later on of all twist-two operators up to two loops 
by means of a conformal operator product expansion (OPE) [26]. 

Here we propose to look at such correlators from a novel point of view. Consider the corre- 
lation function of n protected operators 

G n = {0(x 1 )0(x 2 )...0(x n )). (1.3) 

As long as we maintain the points X{ (with % — 1, . . . , n) in generic positions, this function is well 
defined and has conformal symmetry. As a consequence, it is given by a product of free scalar 
propagators times some (coupling dependent) function of conformal cross-ratios x 2 ^x\ x j \x\x\x) ■ 
In perturbation theory this function is expressed in terms of conformally invariant space-time 
loop integrals. Now, imagine that we wish to take the limit in which the neighboring points 
become light-like separated, 2 

The correlator G n becomes singular in this limit. The first problem we have to face are the pole 
singularities in G n , due to the propagators l/x 2 i+1 connecting two neighboring scalars. Secondly, 
the loop integrals develop logarithmic light-cone divergences ~ lna; 2 j +1 when the integration 
points approach one of the light-like segments [a^Xj+i]. To deal with the first problem, it is 
sufficient to consider the ratio G n /G n \ in which the pole singularities are removed. The second 
problem is more serious, it requires introducing an appropriate regularization. 

Two possible choices of a regularization procedure are: (i) use the small distances xj i+1 as a 
cutoff; (ii) employ standard dimensional regularization and set xf i+1 = from the very beginning. 
These two regularizations are considered in the parallel publication [28], where it is shown that 
in both cases the correlation function reduces to a Wilson loop, 

lim GJGf oc (W[Cn]f ■ (1.5) 

The exact form of the proportionality factor in the right-hand side of this relation depends on 
the regularization; for case (ii) it is just 1. Here W^[C n ] is the light-like polygonal Wilson loop 
in the fundamental representation of the gauge group described earlier. Since H^[C„] is dual to 
the MHV gluon amplitude, Eq. (1.1), we expect that the ratio of the correlation functions in the 
left-hand side of (1.5) is also related to the ratio of amplitudes, A n /A n . The question arises if 
we could find another, more direct way of establishing the relation between correlation functions 
and amplitudes without invoking Wilson loops. 

In the present paper we propose a scenario which realizes this direct link. It employs an 
unusual, dual infrared dimensional regularization procedure. This may seem surprising, as we 
have just argued that the singularities of the correlator occur at short distances. Nevertheless, 
we can consider the following alternative. We start by computing the loop corrections to the 
correlator by means of Lagrangian insertions. This method is well known in field theory and has 
been further developed in [29, 23] for the two-loop computations of four-point correlators, but 
it has universal applicability. The idea is to interpret the loop corrections to the correlator as 
derivatives with respect to the coupling g. For instance, the one-loop correction 

r) C 

9 2 jf^ G n = ~i / d D x Gnl 1 (x ;xi,...,x n ) (1.6) 



2 A similar light-cone limit has extensively been studied in QCD, see, e.g., the review [27]. 



is calculated from the Born-level (n + l)-point correlator 

G^li (x ;x 1 ,...,x n ) = (L(x )O(xi) . . . 0(x n )) (1.7) 

obtained by inserting the Lagrangian at the extra point xq. The crucial point here is that the 
correlator (1.7) stays well defined (after dividing it by the tree approximation Gh ) in D = 4 
dimensions, even if we put the outer points x,i on the light cone, but keeping the insertion point 
xq in a generic position. The logarithmic singularities originate from the integration over the 
insertion point in (1.6). Then we propose to regularize this integral by choosing a measure in 
D = 4 — 2e dimensions, with e < 0. This unusual regularization is motivated by the analogy 
between the space-time loop integrals appearing in (1.6), and the momentum loop integrals in 
the gluon MHV scattering amplitude A n discussed earlier. The analogy becomes possible after 
the T-duality transformation (or change of variables in the integrals) (1.2), provided we use the 
infrared-like regulator above. Still, nothing guarantees at this stage that we will find a result, 
not only similar, but identical with a scattering amplitude. Yet, rather surprisingly, this is what 
happens. 

In the present paper we show a number of examples of this new phenomenon. These include 
all the n— points correlators at one loop, and the four- and five-point correlator up to two loops. 
At present we have no explanation why this is so, but if this new duality correlators/amplitudes 
is confirmed, it will provide the natural explanation of the mysterious dual conformal symmetry 
[13] of the loop momentum integrals in all available gluon amplitude calculations (n gluons at 
one loop [30], four gluons up to five loops [14, 15, 16, 17], five [31] and six [8] gluons up to two 
loops). 

The paper is organized as follows. In Sect. 2 we show a simple example of a correlation 
function becoming an amplitude in the light-cone limit and we formulate the main idea of the 
relation correlators/amplitudes. In Sect. 3 we give a more detailed description of the correlators 
of protected half-BPS operators in superspace. We explain how their loop corrections can be 
obtained by Lagrangian insertions. We then outline the procedure of establishing the duality 
with amplitudes, in particular the introduction of dual infrared regularization. In Sect. 4 we 
show in detail how the duality works for n— point correlators of bilinear half-BPS operators at 
one loop. In Sect. 5 this is expended to four-point correlators up to two loops, and in Sect. 6 to 
five-point correlators up to two loops. Appendix A summarizes some key points of the harmonic 
superspace formalism which we employ for loop calculations. In Appendix B we generalize our 
findings to four-point correlators of half-BPS operators of arbitrary weight. 

2 A simple example of the duality correlators/amplitudes 

In this section we discuss some general properties of the protected operators in M = 4 SYM and 
of their correlation functions. To illustrate our main idea, we show how the four-point one-loop 
correlator is transformed into the four-gluon one-loop amplitude in the light-cone limit. 

The gauge invariant composite operators in M = 4 SYM can be classified as representations 
of the super conformal group PSU(2,2/4) (see, e.g., [32] and references therein). There are two 
basic types of such operators, usually referred to as protected (or "short") and unprotected 
(or "long"). The former satisfy conditions of BPS shortening, i.e. they are annihilated by a 
fraction of the supercharges. This, together with the conditions that they are super conformal 
primaries, implies that they have quantized, or protected conformal dimension equal to their 



canonical dimension. The long operators correspond to generic superconformal representations, 
they receive quantum corrections and acquire anomalous dimensions. 

In this paper we will consider only operators of the half-BPS type. Their lowest components 
(or superconformal primaries) are made of the six real scalars in the M = 4 vector multiplet, 
4>ab = —<Pba = \^abcd4> cd , where A,B = 1, . . . , 4 are indices of the fundamental irrep of the 
R symmetry group S77(4). Generically they are of the type O^ = Tr(0 fc ), carry S77(4) Dynkin 
labels [0, k, 0] and have fixed conformal dimension d = k. The best known example is the simplest, 
bilinear (k = 2) operator, belonging to the so-called stress-tensor superconformal multiplet. The 
top spin state in this multiplet is the stress tensor, while the state of highest dimension is the 
Lagrangian of the M = 4 SYM theory. The lowest dimension state of the multiplet is the bilinear 
scalar operator 

Oabcd = Tr(0 AB CD ) - —e A BCD^{4> EF (t>EF) (2.1) 

belonging to the irrep 20' = [020] of 577(4). Here (J>ab = ^ab^i where t a are the generators of 
the fundamental representation of the gauge group SU(N C ), normalized as tr(t a t b ) = ^5 ab . In 
what follows we always assume the planar limit, 

a = THT , N c -> oo . (2.2) 

Let us consider certain projections of (2.1), namely 

O = Tr(0 12 12 ) , 6 = Tr(0 12 12 ) , 6 = 2 Tr(0 12 12 ) - ^Tr(0 £F £F ) , (2.3) 

o 

where O is the (complex) highest-weight state, O is the conjugate lowest-weight state and O is 
a real projection. We want to evaluate the correlator of n such operators. For n = 2m we can 
take, e.g., m operators O and m conjugates O and consider the correlator 

G n = (0(x 1 )6(x 2 ) . . . 0(x n _i)6(x n )) • (2.4) 

For n = 2m + 1 we can add one operator O, replacing (2.4) by 

G n = (0(x 1 )6(x 2 ) . . . 0(x n _ 2 )6{x n ^)6{x n )) . (2.5) 

Such correlators are finite (the operators (2.3) are not renormalized) and conformally covari- 
ant, as long as the points Xi are kept apart, x; t ^ Xj. If we let two points get close to each other, 
Xi — > Xj, we are dealing with the well-known short distance expansion of the product of operators 
0(xi)0(xj) mentioned above. Here we plan to do something else. We wish to take the limit 
where the neighboring points become light-like separated, without coinciding with each other 3 

x\ i+l ->■ , Xi ^ x i+l , (i = l,...,n), (2.6) 

(with the cyclic condition x n+ i = Xi). This limit is singular for two reasons. Firstly, the 
correlator develops pole singularities, as can be seen already from the (connected, planar) tree- 
level approximation 

Gf = ^-^ c + subleading terms . (2.7) 

x 12 x 23 • • • x nl 




>X n -\ %2\ \ }%n-l 



X n X\ X n 

(a) (b) 

Figure 1: Feynman diagrams of different types contributing to the correlator (2.4) at tree level. Arrowed 
lines denote free scalar propagators {4> l2 {xi)(f)i2{xj)). In the light-cone limit xf i+1 — > the leading con- 
tribution comes from diagram (a), while that of diagram (b) is suppressed by the factor x 2 i x 2 n /(x 3n x 2 4 ). 



By "subleading" we mean terms corresponding to different Wick contractions of the scalar fields 
4> which are less singular in the limit (2.6) (see an illustration in Fig. 1). This can be remedied 
by considering the ratio 

lim GJGSK (2.8) 

Notice that the limit (2.6) breaks the symmetry of G n under the exchange of identical operators 
(e.g., exchanging all points with odd or with even numbers in (2.4)). Instead, it has a cyclic 
symmetry, Xi -> x i+ i, and a flip symmetry, X; L — > x n -i+i. 

Secondly, the loop integrals develop additional light-cone singularities in the limit (2.6). To 
illustrate this, let us take a simple example - the four-point correlator (2.4) at one loop. It has 
been computed in [20, 21, 22] and the result for the ratio G^jG\ is given by 

G A /Gf = 1 + 2a x 2 13 x 2 24 g(xi) + 0(a 2 ) . (2.9) 

Here the one- loop integral g(xj) is defined by 

d 4 x 
z/i j x 10 x 20 x 30 x 40 



9(%i) = — / ~2 ~2 J „2 • ( 2 - 10 ) 



As long as the outer points are kept in generic positions, x 2 i+1 ^ 0, this integral is finite and 
conformally covariant in four dimensions. This allows us to write it down as a function of two 
conformal cross-ratios 

9(x i ) = ^ r *V(%v), u = X ^ff, v = ^, (2.11) 

X 13 X 24: x 13 x 24 x 13 x 24 

where the two- variable function $' 1 ' can be found in [33]. The cross-ratios vanish in the limit 
(2.6) and this function develops a logarithmic singularity, ^'(u, v) ~ lnt> ln-u as u,v — > 0. It 
originates from the integration in (2.10), when xq approaches one of the four light-cone segments 



3 The standard OPE is done in the Euclidean regime, where xf i+1 — implies xi — Xi+\. The Minkowski 



regime allows us to consider the new possibility (2.6). 



Therefore, to define the integral (2.10) in the limit (2.6) we have to introduce a regularization. 
The standard approach is to use dimensional regularization from the very beginning, that is, to 
repeat the whole calculation that leads to (2.9), but in D = 4 — 2e dimensions (with e > 0). This 
approach was adopted in [28], where it was shown that the limit (2.8) turns the correlator into 
a light-like Wilson loop. 

Alternatively, we might use the four- dimensional result (2.9) and declare that we simply 
regularize the integral (2.10) by choosing a measure in D = 4 — 2e dimensions (with e < 0). 
Notice the change of sign of the regulator - now it looks more like an infrared, rather than the 
natural ultraviolet regulator needed for such short-distance singularities. This unusual choice 
is motivated by the observation that the one- loop space-time integral (2.10) is the dual space 
version [13] of the one-loop scalar box momentum integral appearing in the four-gluon MHV 
amplitude. The latter is given by [15] 

A 4 /4 0) = 1 + astll 1 \p i ) + 0(a 2 ) , (2.12) 

where A\ is the tree-level amplitude, s = (pi +P2) 2 , t = (ps + Pa) 2 are the Mandelstam variables 
and 

/«(„.) = 2 * f ft (2 13) 

4 {Pi) (27iy-^Jk^k- Pl y(k-p 1 -p 2 y(k + p 4 y y ' } 

is the one-loop scalar box integral. Switching from momenta to dual coordinates, k = X\ — Xq 
and pi = Xi — Xi+i (with x 5 = xi) , we identify the two integrals: 

li%i)=9 e (xi), (2-14) 

where the subscript in g e indicates that we have changed the measure in (2.10) to D = 4 — 2e 
dimensions (with e < 0). Notice that the light-cone limit (2.8) for the correlator implies that the 
momenta pi = x^i+i are light-like, p 2 = 0, as required for the amplitude. 

Further, comparing (2.9) and (2.12), we observe a surprisingly simple relation between cor- 
relator and amplitude: 

lim G 4 /Gf ] = (A 4 /Af ] ) 2 + 0(a 2 ). (2.15) 

This one-loop exercise suggests the following general recipe for obtaining MHV gluon ampli- 
tudes from correlators of protected operators: 

1. Compute the n— point correlator of protected operators in four dimensions. 

2. Change the D = 4 integration measure in the loop integrals to D = 4 — 2e (with e < 0). 

3. Divide the correlator by its tree-level value and take the light-cone limit x 2 i+1 —$■ 0. 

4. Switch from coordinates to dual momenta. The result should be the square of the n— gluon 
MHV amplitude, divided by the tree. 

Of course, the four-points/gluons case at one loop is probably too simple to allow us to jump 
to the conclusion that such a duality correlator/amplitude is a general phenomenon. It is the 
purpose of this paper to provide a lot of evidence in favor of this conjecture. 



3 Correlators of protected operators in superspace and 
their loop corrections 

In this section we explain the role of the Lagrangian insertion procedure in the calculation of loop 
corrections to correlators. Then we formulate the rules of dual infrared dimensional regularization 
and state the main result of the paper. 

3.1 Correlators of protected operators in J\f = 2 harmonic superspace 

The one-loop, and even more so the two-loop calculation of the correlators are much easier to 
do in terms of superspace Feynman graphs. Since we do not have an off-shell formulation of the 
M = 4 SYM theory suitable for application to perturbation theory, the best compromise is to 
use the formulation in terms of M = 2 superfields in the so-called harmonic superspace [34, 35]. 
We give a brief summary of this formalism in Appendix A. Here we just mention the types of 
M = 2 supermultiplets and superfields that we are dealing with. 

The M = 4 vector multiplet is decomposed into an M = 2 matter multiplet (hypermultiplet) 
and an J\f = 2 vector (gauge) multiplet. Upon reducing the R symmetry, SU(4) — > SU{2) x £7(1), 
the six real scalars <Pab split into an isodoublet 0* (with % — 1,2) and a complex singlet </?; the 
four (chiral) gluinos \ Aa split into a doublet \ ta and two singlets ip a , K a (and their antichiral 
conjugates). These fields can be combined to form M = 2 superfields. One of them describes the 
hypermultiplet: 

q + (x,9 + ,9 + ,u) = ( p i (x)u+ + 6 +a ^ a (x) + 9+R & (x) + ... (3.1) 

(the dots denote auxiliary and derivative terms). It is a Grassmann analytic (or G-analytic, or 
half-BPS) superfield in the sense that it depends on half of the Grassmann variables, 9 +a = 9 ta uf 
and 9 +a = 9 ia uf, obtained by projecting the SU(2) doublets 9 ta and 9 ta with an SU(2) harmonic 
variable uf. The latter, together with its conjugate u~ = (u +1 )* forms an SU(2) matrix. Notice 
the presence of both chiral and antichiral odd variables in the hypermultiplet superfield (3.1). In 
harmonic superspace one can define a special conjugation (denoted by tilde ~), which takes (3.1) 
to another G-analytic superfield, 

q + (x,9 + ,9+,u) = 4>\x)ut + +a K a {x) + 9^(x) + . . . . (3.2) 

In contrast with the hypermultiplet, the M = 2 vector multiplet is described by the chiral field 
strength (and its antichiral conjugate) 

W(x, 9) = (p(x) + 9 ia \ ia {x) + 9 ia 9?F a p(x) + ... (3.3) 

containing, in particular, the self-dual part of the gluon field strength F a/3 = ( y cr^ u ) al3 F fJiU . 

The protected half-BPS bilinear operators (2.1) in the SU(4) irrep 20' split into a number 
of irreps of SU(2) x £7(1). They can be descried as the lowest (9 = 0) components of bilinears 
made of the above superfields. Three such bilinears are G-analytic superfields: 

= Tr(q + q + ), = Tr(q + q + ), 6 = 6 = 2 Tr(q + q + ) , (3.4) 

where all operators are functions of x, 9 + , 9 + , u. For example, the operator 0(x, 9 + , 9 + , u) has the 
bottom component 0(x,u) = O\ e + = §+ =0 = Tr (0* (x) (j>> (x))ufu^ containing the complex SU(2) 



triplet Tr(0*0 J ). Another, real triplet is the bottom component of the real operator = 0, 
0(x,u) = O\ e+= §+ =0 = 2Ti(4> t (x)(f) : '(x))u'lu^. We may say that the harmonic variables serve as 
a "bookkeeping device" for organizing the fields into £77(2) representations. 

The remaining SU(2) x U(l) projections of the 20' are described by different types of su- 
perfields. Among them, the chiral operator Tr(WW) plays the prominent role of the M = 2 
SYM Lagrangian. In this paper we do not consider the rest of the M = 2 projections of the 20': 
Tr(WW), Tr(WV), Tr(Wq + ) and conjugates. 

The correlators of the three G-analytic operators (3.4) will be the main subject of this and 
the following sections. We are considering n-point correlators of the type, e.g., 

G n = (0(x 1 ,u 1 )0(x 2 , u 2 ) • ■ • 0(x n -i,u n -i)0(x n , u n )) (3.5) 

for n even, or 

G n = (0(xi,ui)0(x 2 , u 2 )--- 0(x n - 2 , u n - 2 )6(x n -i, u n - 1 )6(x n , u n )) (3.6) 

for n odd. At tree level the connected contribution to the correlator has the form 4 



Nl (12)(23)---(n l) 

('0 Tr "\2n rpl rf.2 „2 



G(°) = - c K ! K I? + subleading terms , (3.7) 



where (r, r + 1) is a shorthand for the SU(2) invariant but U(l) covariant contraction of the two 
harmonics with labels r and r + 1. 



( r> r + 1) = -(r + 1, r) = u^jU^ . (3.8) 

In (3.7), as in the expression (2.7) for the tree amplitude for the highest-weight and lowest- 
weight state projections of the 20', we only show the leading singular term in the light-cone limit 
(2.6). Compared to (2.7), the tree-level expression (3.7) contains additional information about 
the isotopic £77(2) structure of the correlator. It is carried by the harmonic variables u^ % (with 
% — 1, 2 and r = 1, . . . , n) at each point. 

3.2 Lagrangian insertions 

The correlators (3.5) and (3.6) are defined by the path integral 

G n = /"p$ e ^= 4SYM 0(x 1 ,u 1 )...0(x n ,u n ) 

= G^ + g 2 GU + g * G ^ + . . . . (3.9) 

Here the M = 4 SYM action consists of two parts, the M = 2 SYM action and the action of the 
Af = 2 hypermultiplet matter coupled to the gauge sector: 

£Af=4 SYM = £aT=2 SYM + £V=2 matter • (3.10) 



4 More precisely, an n-point correlator may involve m complex operators O and also m conjugate operators O, 
the remaining n — 2m operators being of the real type O. At tree level, such a correlator equals c nm G„ , where 
Cno = 1 + (— 1)™ and c nm = (_i)™+ m : m > 0. In this article we are interested in the ratio of loop corrections 
over the corresponding tree, which is universal. The coefficients c nm can thus safely be omitted. 



Instead of computing the loop corrections to G n directly, we prefer to evaluate the derivative 
with respect to the coupling. As we show below, it is given by the insertions of the M = 2 SYM 
action 

SV=2 sym = / d D xd 4 6L Af=2 sym(x,6) = / d D xd 4 8Lj^ =2 sym(x,9) , (3.11) 

where, after the appropriate rescaling of all the fields of the M = 2 vector multiplet, 
1 



N=2 sym = — Tr(^ 2 ) (3.12) 

1 



2g 2 



-^-{Tr(y? 2 ) + . . . - {e) A Ti[F a0 F ap + 0n<f>i + UX'dXi + interaction terms]} 
2g 2 

is the M = 2 SYM chiral Lagrangian. The effect of the rescaling is that the coupling g disappears 
inside the Lagrangian Lj^ =2 sym, it is only present in front of it, as indicated in (3.12). The cou- 
pling also drops out from the interaction of the vector and matter multiplets (see Appendix A. 3). 

Notice that the action in (3.11) has two forms, one chiral, the other antichiral. They are 
equivalent due to a Bianchi identity stating that the difference between Tr(W^ 2 ) and Tr(W 2 ) is a 
total (super) space derivative. At the component level, this is clearly seen from (3.12), where, for 
instance, the complex combination F a pF al3 = F 2 + iFF contains the Yang-Mills Lagrangian F 2 
and the topological term iFF . We will come back to this important point later on. 

Now, we want to differentiate the correlator with respect to the coupling, which is present 
only in SV=2 sym as an overall factor, see (3.12). Thus, the one- loop (order g 2 ) correction to the 
correlator, 

g 2 Tr^ G n = -i d D x Gn+iixo; x ± , wi; . . . ; x n , u n ) + 0(g 4 ) (3.13) 

is calculated from the Born-level (n + l)-point correlator 

Gnl 1 (x ',x 1 ,ui',...',x n ,u n )= / d A 6 (Ltf =2 sym{xo,0 )0{x 1 ,u 1 ) . . .0(x n ,u n )) + 0{g 4 ), (3.14) 



obtained by integrating the Lagrangian insertion over the Grassmann variables at the insertion 
point, but not over the space-time point Xq. Note that this tree-level correlator is of order 0(g 2 ), 
because it involves interaction vertices (see Sect. 4 for details). 

A very important property of the correlator <3 n +i i s i^ s superconformal symmetry. Indeed, 
it involves the protected operators O and L/v"=2 sym (the latter belongs to the N = 4 stress- 
tensor multiplet), with fixed conformal dimensions 2 and 4, respectively. Such operators are not 
renormalized and have well-defined conformal properties. This symmetry greatly facilitates the 
perturbative calculation, as explained in Appendix A. 

The same procedure can be applied to the higher-order perturbative corrections. Thus, to 
obtain the correlator at two loops (order g 4 ), we compute the derivative 

o gA ( Jp2 ) Gn = ~2 dDx ° dDx °' GnU&o, x >] xi, wi; . . . ; x n , u n ) + 0(g 6 ) (3.15) 

in terms of the Born-level (n + 2)-point correlator with two Lagrangian insertions (see Ap- 
pendix A. 3) 

Q^l 2 (x , x r,X!, Mi; ... ; x n , u n ) = / d 4 9 d\,(L(x , Q )L(x >, 9 >)O(x 1 ,u 1 ) . . . 0(x n , u n )) + 0(g & ) . 

(3.16) 



In conclusion, the Lagrangian insertion procedure reduces the calculation of loop corrections 
to the correlator G n to a tree-level calculation of the correlator with insertions. The main point 
of our conjecture is that this tree-level correlator tells us what the integrand of the dual MHV 
amplitude should look like. The precise matching of the two objects is obtain by introducing a 
dual infrared regulator. 

3.3 Outline of the dual infrared regular izat ion procedure and of the 
duality correlator/amplitude 

Our strategy for establishing the relationship between correlators and amplitudes is as follows. 
We start with the tree-level correlators Gn+i (3.14) and Q n l 2 (3.16). They are computed in 

D = 4 and need no regularization. Then, we put Gn+i an d G n +2 on the light cone by setting the 
adjacent external points at light-like distances, xf i+1 — > 0, i = 1, . . . ,n. In order to remove the 
pole singularities, we divide the correlator with insertions by the tree-level correlator without 
insertions (3.7), thus obtaining U(l) chargeless ratios, Gnl\/G n and Q n \. 2 IGn ■ After that it 
becomes safe to set the external points on the light cone, while keeping the insertion points xq 
and xqi in arbitrary positions. We remark that at this stage we still need no regularization. 

Next, we perform the integration over the insertion points, thus passing from the tree-level 
correlators with insertions Gn+i an d G n +2 t° the loop corrections of the correlator without inser- 
tions G n . Here we are facing logarithmic singularities due to the divergent integrals over xq and 
xqi. Divergences arise when the integration points approache a light-like segment [xj,Xj + i]. We 
regularize the integrals by modifying the dimension of the integration measure, D = 4 — 2e, with 
e < 0. We emphasize that this is not standard dimensional regularization, for two reasons: (i) 
the tree-level correlator with extra points from the Lagrangian insertions has been computed in 
D = 4 and then put on the light cone, without regularization; (ii) the sign of the regulator e is 
chosen to match the infrared divergences of the dual amplitude, so this is not the usual ultraviolet 
regulator. We call this "dual infrared regularization" . 

In order to make contact between the n-point correlator and the n-gluon amplitude, we 
identify the momenta with the dual coordinates 

Pi = 2%i+l j x n+ \ = X\ , (3-17) 

so that YliPi = an d Pi — 0. The Mandelstam variables are identified with the non- vanishing 
distances in dual space, 

Sij = (Pi + Pi+i + ■■■+ Pjf = x\ j+1 . (3.18) 

The main result of the next three sections is that the above procedure leads to the following 
duality relation between correlators restricted to the light cone and MHV n-gluon scattering 
amplitudes: 

lim In {GJGf) = In {AJA^f + 0(e) , (3.19) 

where A n is the tree- level amplitude and 0(e) denotes terms that vanish after we remove the 
regularization. The reason why we formulate the relation in terms of logs will become clear in 
Sect. 6. In Sect. 4 we demonstrate this duality for any n at one loop, and in Sections 5 and 6 we 
will show it for n = 4, 5 up to two loops. 
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4 From correlators to amplitudes: n points at one loop 

In this section we perform the calculation of the one-loop correction to the correlator G n ac- 
cording to the procedure of Sect. 3.3. 5 We use the Feynman rules from Appendix A. 2. The 
harmonic superspace Feynman diagrams for the correlator with one extra point corresponding 
to the Lagrangian insertion are shown in Fig. 2. 6 



fc + 1 





(a) 



(b) 



Figure 2: The n-point correlator with one insertion. The solid and wavy lines are hypermultiplet and 
gauge propagators, respectively. The double dot denotes the insertion of the N = 2 SYM Lagrangian 
Tr(FT 2 ). 

They are constructed from the two basic building blocks T and TT described in Appendix A. 3. 2, 
combined with free hypermultiplet propagators (A. 24). The block T shown in Fig. 3(a) is the 
supersymmetric analog of the vertex correction {4> a {xi)F{^(xo)(j) c (x2)) for two scalars and one 
gauge field strength. The main difference is that at the insertion point we have only the self- 
dual part of the gauge field strength F a p = {cr^apF^, as well as the auxiliary field Y^ (see 
(A. 13)); the scalar (p and the gluino X ta from (3.3) do not have a cubic vertex with the external 
scalars. Since we are doing the calculation in D = 4, the integral at the vertex can be easily 
computed yielding a rational expression. It is obtained from the block (A. 47) by setting the 



external 9^ = 9 2 



0: 



(?+(xi, 0, Ui)W b (x , 9 )q+(x 2 , 0, u 2 )) 



2ig 2 f abc 



12) 



(2tt)< 



X 



-%\2 



(4.1) 



12 



with 



7-12 



X 



12 



7 0/l ' u 0/2 

12)x 10 x 2 q 



y 0/l ' u 0/l 



'0/2 ' u 0/2 



X 



10 



X 



2(1 



7 /lF 10 1^20^0/2 

\Y2)x w x 2 q 



(4.2) 



Here 9^, = 9 l (u r )f are the two U(l) projections of the SU(2) doublet 9 % with the harmonics at 
point r; 9 ■ 9 = 9 a 9 a ; [x,y]J = x a6l y^ - y a& x*P = -2ix^(a, 



\0 



' iiV) a 



5 Our result agrees with earlier calculations of one-loop correlators in [20, 21, 22, 36]. 

6 Notice that the graphs are drawn with a polygonal matter frame. Graphs based on the "zigzag" configurations 
like in Fig. 1(b) are suppressed in the light-cone limit (2.6), after dividing out the leading singularity of the tree- 
level correlator (2.7). 
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(a) (b) 

Figure 3: Building blocks for the graphs in Fig. 2. The single dot in (a) denotes a W insertion, the 
double dot in (b) a Tr(VF 2 ) insertion. 

The second building block shown in Fig. 3(b) is the supersymmetric analog of the propagator 
correction ((f>(xk)<f>(xk+i)) > after integration over the insertion point. It is obtained from the TT 
block (A. 48) by identifying the insertion points = 0' and by setting 9\ = 9\ = ; 

(qt(x 1 ,0,u 1 )Tr(W 2 )(xo,e )qt(x 2 ,0, U2 )) = 2 ffi** flftj 12 (4.3) 

{Z7T) X 12 

with 



; i-2-) W n ?' a + 2 

(12) ^10^20 



J 12 = 777^ 2 2 _ 2^ 12 ) ■ ^ 4 ' 4 J 



Summing up all graphs (all sums are cyclic, k + n = k) and integrating over the odd variable 
#o at the insertion point (but not yet over the even xq), we find 



n n 



c 2 _IK (12) (23) •••(»!) I" ,„ (1 ^ . . ^.. n 

^+! - / 97r ^2(n+2) T 2 T 2 3 / " *° U 2^ *fc,fc+l*I,*+l + 2^ ■?*. fc + 1 J 

^tt; x 12 x 23 x nl j KZ ktl=1 . ik¥i fc= i 7 

= t/^1°) jd%^ Kk+1 )\ (4.5) 

71 ^ fc=i 

where G*i is the leading singular part of the connected tree-level n-point correlator defined in 
(3.7). In deriving (4.5) we have used the identity ju — \i}ki) 2 ■ 

The next step is to set the adjacent external points on the light cone, x\ k+1 = 0, for the ratio 



/n 
d% (J2i k , k+ 



i, 2 . (4.6) 



We notice that for x\ 2 —* the first term in i\2 (4.2) vanishes while the second and the third 
terms cancel in the cyclic sum 



fc=i 
Then, using the identity 



EA^ /fc ko,%i#J fc+1 
^ fc +! - ~ 1^ (u h + l)r 2 r 2 ' ^ ' 

k=l \r», ft- t 1 Mfc,o x fc+i,o 



/ d% 9% lh 9% 2i2 9% 3h 9% 4i4 = -- ( e i l i 2 e '3*4 e a 1 a 4e a 2 a 3 _ ^4^3^1*2 ^a^ ^ ^ 
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after some algebra we obtain 

2 

lim 



X 2 



^ oJ XT J 

= 2 lim n L~J^ ir~3 / d4 ^o^J fc [aJfc,0' a: fc+i,o]^/fc + i^J'j[a ; J,o,aJj+i,o]^Ji + i 



fc j **'fc,0" i/ fc+l,0**'i,0**'l+l,0 



lim 



{^ 



( x k,0 ■ Xl+l,0)(Xk+l,0 ' x lfl) — { x kfl ■ x l,o)( x k+l,0 ■ Xl+lfi) 



x tt+i-^° *• fc j • E fc,0 a 'fc+l,0 a 'J,0 a 'J+l,0 



+'E 



A:./ 



t 1 v X P 

t ^^P- L k,0- L k+l,0- L lfl- L l+l,0 i / . „s 

x fe,0 x fc+l,0 x «,0 x i+l,0 



The parity-odd terms in the last line deserve a special comment. They will subsequently be 
integrated over the insertion point xq, according to (3.13). The resulting pseudo-scalar integral 
will depend on the four external points Xk, Xk+i, xi, xj+i, which is not sufficient to make a transla- 
tion invariant pseudo-scalar. But in fact there is another reason why these terms are suppressed 
by the integral over the insertion point. By inspecting the component content of the inserted 
Lagrangian (3.12), we see that the pseudo-scalar terms are due to the presence of the topological 
term IFF. Such an insertion is a total space-time derivative with respect to the insertion point, 
so the integral must vanish. The underlying reason why the loop corrections to the correlator 
G n cannot contain parity odd terms is that the fields in the operators O can be treated as true 
scalars (see Appendix A. 5). We will come back to this point in Sect. 6. 

Further, replacing 2(xk,o • %i+i,o) — x \ o + x f+i o~ x t i+n e ^c in (4.9) and using the properties 
of the cyclic sum over k and I, we obtain 



2 



A^a fV- \ L v V^ x k,l+l x k+l,l x kl x k+l,l+l (Ain\ 

d 9 2^i k>k+1 = - km l > -5-5 . (4.10) 

\ fc / Zx t>+^°k~l X k,0 X k+l,0 X l,0 X 1+1,0 

Notice that the numerator vanishes when k — I — 1, I or I + 1. 

Up to now, we have done the entire tree-level calculation, including the light-cone limit, in 
D = 4. The last step is to integrate over the insertion point, J d D xo, with D = 4 — 2e and e < 0, 
in order to regularize the divergences occurring when the integration point approaches any of 
the light-like segments [xj,Xj+i]. We reiterate that this is not the natural UV regularization (for 
which e > 0), but it is what we call a dual infrared regularization. This yields 

r & (r< /ri(0)\ i i- V^ f jD x k,l+l x k+l,l ~ x kl x k+l,l+l //iii\ 

2 lim _n ITa (° n/Gn > = ~i^ 2 hm >ni- / d X ° r 2 r 2 r^T 2 ' (4U) 

x i,i+i-+° ° a a=0 47r X U+^°T X k,0 X k+l,0 X l,0 X l+l,0 

To make a comparison with the amplitudes, it is instructive to rewrite the one-loop integral 
in the right-hand side of (4.11) in terms of dual momenta, Eqs. (3.17) and (3.18). Defining the 
D— dimensional loop momenta as £ = xo,k we find that, for general k and I this space-time integral 
becomes the two-mass easy box momentum integral [30] 

F(p P g Q) - --L f fi W ~ (P + P ^ + ^) 2 ) (p 2 = , 2 = 0) (4 12) 
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evaluated for p = Xk,k+i, P — %k+i,i, Q = %i,i+i and Q = xi+i t k- When rewritten in terms of 
these functions, the right-hand side of (4.11) coincides with twice the one-loop n-gluon MHV 
amplitude (see Eq. (4.19) in [30]) 

.Mm |-(<?„/Gf) =24 1 )/4°) = 2^F(p,P,g,Q). (4.13) 

In conclusion, the correlator calculated by means of a Lagrangian insertion and in dual infrared 
dimensional regularization, reproduces the amplitude in terms of momentum integrals. 

5 Four-point correlators and four-gluon amplitudes 
to two loops 

In this section we extend the one-loop duality from the previous section to two loops. We show 
that the logarithm of the correlator of four half-BPS operators of weight two (bilinears), when 
put on the light cone using the dual IR regularization procedure, becomes identical with the 
logarithm of the square of the four-gluon scattering amplitude. 

5.1 Four-point correlators of bilinear half-BPS operators 

Let us consider the correlator of four protected M = 2 half-BPS complex operators (3.4) 

G A = {0(xi,ui)0(x 2 , u 2 )0(x 3 , u 3 )0(x 4 , u A )) . (5.1) 

It has been computed up to two loops in [23] 7 (we summarize the computation in Appendix A): 8 



3 ( 0) = N 2 (12) (23) (34) (41) 

( / ijr \ O T iy>£ rp£ iy^ 

\4n) X 12 J-23 X 34 X 



^4 _ /o_\S _2 _2 ..2 ..2 ' \°- Z ) 



G 4 — \ s — — 2 — 2 — T x i3 a; 245 , o(l,2,3,4) , (5.3) 



(1) = N 2 2aU , 2 

( O/TTlO 71* /-**£ /y*2> /y*£ 
\&ii ) X-L2X23X34X4 

( 2 ) = N 2 2a 2 TZ 

yz,li) J -12 X 23 X 34 X 41 



yz,H) J-i2 x 23 J '34 J '41 

— (£12^34 + X l , i X 2A ' X 14 X 23)\90{1-, A 3, 4JJ 



+ x 2 13 x 2 M (x 2 12 h (l, 2, 3; 1, 2, 4) + x 2 23 h (l, 2, 3; 2, 3, 4) + x 2 34 h (h 3, 4; 2, 3, 4) 
+ x 2 41 h (h 2, 4; 1, 3, 4) + x 2 13 h {^ 2 , 3; 1, 3, 4) + 4^(1, 2, 4; 2, 3, 4) 



(5.4) 



7 The terminology of [23] is based on the topology of the Feynman graphs, rather than on the perturbative order 
of the correlator. Thus, the lowest order (Born level) contribution is called "one-loop" , and the contributions of 
order g 2 ^ -1 ) arc called "£-loop". 

8 Thc Born approximation (5.2) has a connected and a disconnected sectors (see [23]), with different color 
factors. Here we show only the connected contribution to the tree-level correlator. Also, the color factors are 
given in the large N c approximation. 
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where the one- and two-loop four-dimensional integrals in coordinate space are defined by 

/d 4 x 
T 2 T 2 1 T 2 > ( 5 - 5 ) 

x 10 x 20 x 30 x 40 

/io(l,2,3;l,2,4) = cg / 2 2 2 2 2 , ,,, (5.6) 

J ^10 x 20 x 30M00'v x 10' x 20' J '40 / / 

with Co = i/(2ir 2 ). A characteristic feature of the all-order loop corrections is the presence of the 
universal harmonic-space-time polynomial prefactor 

U = (12) 2 (34) 2 x? 4 a4 + (14) 2 (23) 2 :r 2 2 x 2 4 + (12) (23) (34) (41) [x 2 3 x 2 4 - x\ 2 x 2 u - x 2 4 x 2 3 ] . (5.7) 

This phenomenon was first revealed in [37] under the name "partial non-renormalization" . 

Note that if the external points are in generic positions with a; 2 7^ 0, the integrals above are 
well defined in D = 4 and are manifestly conformally covariant. This allows us to write them as 
functions of the conformal cross-ratios, 



So(l,2,3,4) = ^$W(^,^f ) 

x 13 x 24 \ x 13 x 24 x 13 x 24/ 



Ml, 2, 3; 1, 2, 4) = c 2 / , 2 2 f*\ dD f X " 2 2 , , (5.9) 



-1 2 2 9 2 

Ml,2,3;l,2,4)= *^«P) ¥¥,¥^ , (5-8) 

l x 12/> x 34 vx 12 x 34 x 12 x 34 / 

having well-known expressions in terms of polylogs (see, e.g., [33]). In what follows we wish to take 
the light-cone limit x 2 i+1 — > 0, which makes the integrals diverge. We regularize them, mimicking 
the gluon scattering amplitudes, by modifying the dimension of the integration measure to D — 
4 — 2e (with e < 0), thus giving up conformal invariance 

#(l,2,3,4) = c e / — — - — \ — r , 

J J-io x 20 x 30 x 40 

d D x 

(,y,Z ,y,Z rpA \ ,-pZ I ,y, A ,y, A rpA \ 

x 10 x 20 X 3(V x 00' \ x 10' x 20' x 40'/ 

where the normalization factor c t = 2i/(27r) 2 ~ 2e is introduced to simplify the final expressions 
for the two-loop corrections. 

5.2 Four-gluon amplitudes 

We expect that the four-point correlator (5.1) in the light-cone limit (2.6) is related to the four- 
gluon planar MHV amplitude. The latter has been computed to two loops in [15] and it has the 
following form: 

A 4 = 4 0) [1 + aM {l) + a 2 M {2) + 0(a 3 )] , (5.10) 

where A 4 is the tree-level amplitude and the loop corrections M^ and M^ are given by 

M (1) = stli 1] (s,t) , M (2) = st(sli 2) (s,t) + £/f } (t,s)) , (5.11) 
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where s = (p± + P2) 2 and t = (p2 +P3) 2 are the Mandelstam variables. All gluons are outgoing, 
so that their momenta satisfy the relations Yli Pi = an d P 2 = 0. 

The one- and two- loop scalar box momentum integrals I A and l\ in (5.11) are given by 



li 1 \s,t)=c e J 
lf\s,t)=c 2 



d u k 



k 2 {k — p\) 2 {k — pi — P2) 2 (k + Pi) 2 ' 

d D k d D l 



k 2 {k — pi) 2 (k — p\— P2) 2 ik + l) 2 l 2 {l — Pi) 2 (l — Pz — Pa) 2 



(5.12) 
(5.13) 



These massless integrals are infrared divergent, therefore they are regularized dimensionally, 
D = 4 — 2e with e < 0. Switching from momenta to dual coordinates [13], k — X\ — xq, 
I — xo> — x\ and Pi = Xi — Xi+i (with x§ = x\) , we identify the pseudo-conformal integrals from 
the correlator with those from the amplitude (see Fig. 4): 

li 1] (s,t)= 0(1,2,3,4), /f ) ( S ,t) = / i (l,2,3;l,3,4), if(i, s) = h(l, 2,4; 2, 3,4) . (5.14) 

Then, the two-loop expression for the amplitude (5.10) can be rewritten as 

A 4 /Af ] = 1 + ax 2 13 x 2 24 g(l, 2, 3, 4) + a 2 x 2 13 x 2 4 [x 2 13 h(l, 2, 3; 1, 3, 4) + x\ A h(l, 2, 4; 2, 3, 4)] . (5.15) 
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Figure 4: One- and two-loop pseudo-conformal integrals contributing to the correlator G4, Eq. (5.1), 
and to the amplitude A 4 , Eq. (5.10). The diagrams with solid lines depict Feynman integrals in x— space. 
The diagrams with dashed lines represent the same integral in the dual momentum space. The straight 
labels correspond to the points Xj, the slanted labels correspond to the momenta pi = X{ — Xj+i. 
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5.3 Duality correlator /amplitude 

For the purpose of comparing correlators and amplitudes, we define the ratio of the correlator 
and its Born-level expression: 

G * /G ^ = 1+ (12)(23)(34)(41) (5 ' 16) 



x<^ 2a x{ 3 x 2 24 g{l, 2, 3, 4) + a 



3 '12 a '34 ~r 3 '13 a '24 ' X 14, X 23)\9{^y 2, 3, 4JJ 



+ 2x{ 3 x 2 24 xi 2 h(l, 2, 3; 1, 2, 4) + x 2 23 h(l, 2, 3; 2, 3, 4) + x 2 34 h(l, 3, 4; 2, 3, 4) 



+ x 2 41 h(l, 2, 4; 1, 3, 4) + x 2 13 h(l, 2, 3; 1, 3, 4) + x 2 24 h(l, 2, 4; 2, 3, 4) 



+ 0(a 3 



Next, we wish to evaluate this ratio on the light cone, i.e. with 

2-12 = -^23 = x 34 = x 41 = U . (O.ll) 

From (5.7) we see that the prefactor 7£/[(12)(23)(34)(41)] in (5.16) is reduced to x\ 3 x\ 4 . Further, 
most of the two-loop h— integrals do not contribute to the right-hand side of (5.16) due to 
vanishing kinematic prefactors like x\ 2 . The result is 



M - 1 %.2.2 



lim G4/Q = !+ 2a xf 3 x 2 24 g(l, 2, 3, 4) + a 2 



*m+i-*> 



(^3^(1,2,3,4))' 



+ 2x\ 3 x 2 24 x 2 13 h(l, 2, 3; 1, 3, 4) + x 2 24 h(l, 2, 4; 2, 3, 4) 



+ 0(a 3 ). (5.18) 



Comparing this result with the perturbative expansion of the amplitude, Eq. (5.15), we obtain 
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lim G 4 /G { 4 0) = ( A 4 /A { 4 0) ) +0(a 3 ), (5.19) 

confirming the general statement (3.19). 

We recall that G4 was defined in (5.1) as the correlation function of operators (2.3) bilinear 
in the scalar fields. The duality (5.19) can be extended to correlators of half-BPS operators 
O = Tr(<f) k ) of arbitrary weight k, see Appendix B. 

6 Five-point correlators and five-gluon amplitudes 
to two loops 

In this section we extend the previously found duality between dual-IR-regularized correlators 
and gluon MHV amplitudes to the case of five points/gluons. This is a rather non-trivial test, in 
view of the significantly more complicated integrals involved. Also, we explain why the duality 
should be formulated in terms of logs, rather than the correlator/amplitude themselves. The 
reason is in the parity-odd (pseudo-scalar sector) of the amplitude, which is reduced to 0(e) 
terms by taking the log. This is essential for the duality to work, because the correlator must be 
a true scalar, as we argue in Appendix A. 5. 
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We consider the correlator of n = 5 half-BPS operators of weight k = 2, 

G 5 = ((9(1)0(2)0(3)0(4)0(5)) , (6.1) 

where the local scalar operators O, O and O are the bottom components of the hypermultiplet 
bilinears 

= Tv(q + q + ), = Tr(q + q + ), 6 = 2Ti(q + q + ), (6.2) 

and i = (xi,uf) denotes the set of space-time and harmonic coordinates of the hypermultiplet 
scalar fields, q + (i)\g=o = (p r (xi)uf r and q + (i)\e=o = (p r (xi)uf r (with r — 1, 2). We wish to examine 
the correlator (6.1) in the limit 

X i,i+1 = \%i x i+l) ^ U , x i+5 = x i ■ (,"•"/ 

In this limit, the leading asymptotic behavior of the tree-level correlator is given by the product 
of free scalar propagators 

(0) _ N 2 (12)(23)(34)(45)(51) 

^5 /o \10 „2 2 „2 2 „2 T • • • , lO-^j 

^Z/i; X 12 J '23 J '34 X 45 J '51 

where the ellipses denote subleading terms as 4+i — > 0. Notice that the obvious symmetry of 
the correlator (6.1) under the exchange of operators, 1 f> 3 and 2 -h- 4, is lost in the light-cone 
limit. At loop level, the correlator (6.1) turns out to have the same leading singularity as the 
tree G 5 . This suggests to study the following ratio in the light-cone limit (6.3), 



F G = lim In (G 5 /G£A = aF%\ Xi ) + a 2 F%\ Xi ) + 0( 



a 6 ) , (6.5) 



where a = g 2 N c /(8ir 2 ) and Fq (xi) are scalar functions of Xi only. The rationale for considering 
the log of the ratio of correlators in the left-hand side of (6.5) is that, firstly, it does not receive 
a contribution at 0(a°) and, secondly, as we will argue below, it has a much simpler form. 
Computing Fq(xi) and Fq (xi) we shall follow the same routine as before, that is, we shall 
expand Fq (xi) over a basis of one- and two-loop pseudo-conformal integrals in D = 4 dimensions, 
and then regularize them by modifying the integration measure at the Lagrangian insertion points 
to D = 4 — 2e dimensions, d 4 x — > d 4 ~ 2t x , with e < 0. 

In Sect. 4 we have demonstrated that the one-loop correction in (6.5) is given by 

F G 1 \x i ) = 4^(1,2,3,4) + (cyclic) , (6.6) 

where #(1,2,3,4) is the one-loop pseudo-conformal "cross" integral in (6.8); (cyclic) means the 
four non-trivial cyclic permutations, % \- > i+1, of the points {1, 2, 3, 4, 5}. The two-loop correction 
F G is computed in Appendix A. 4, using the method of double Lagrangian insertions in harmonic 
superspace. The result is the following expression: 

F G 2) =x 2 13 x 2 24 [2 x 2 n h(l, 2, 3; 1, 3, 4) + 2 x\ A h(l, 2, 4; 2, 4, 3) - x 2 u h(l, 2, 4; 1, 3, 4) 
+x 2 13 x 2 u x 2 25 [2h(l,2,5;l,3A) ~ h(l,2,3;l,A,5) - h(l, 2,4; 1, 3, 5)] 



+4^35 [2 4 5 p(l;2,5;3,4) - x 2 24 p(l; 2,4; 3, 5) - x 2 35 p(l; 3, 5; 2,4) 

-\[xl 3 x 2 2i g(l, 2, 3, 4)] 2 - 4 4 4 «,(!, 2, 3, 4) 4 x 2 25 g(l, 2, 3, 5) + (cyclic) , (6.7) 



where the cyclic permutations act on the entire sum of terms. Here we use the dimensionally 
regularized integrals (5.9) together with 

p(l; 2, 3; 4, 5) = c\ \ - 2 ^ 2 f^f^ f - 2 - 2 ; , (6.8) 



(rp^ rp^> rp^ \ rp£ I ,-v* rp^ rp^ iy^> 

x 10 x 20 X 30J x 00' v x 20' x 30' ^W ^50' > 

where D = 4 — 2e with e < and the normalization factor c e = 2i/(27r) 2_2e is introduced to 
simplify the final expressions for the two-loop corrections. Switching to dual momenta pi = 
Xi — Xi+i, we find that the g—, h— and p— integrals correspond to the scalar box, double-box and 
penta-box momentum integrals of Ref. [31], respectively. 

6.1 Five-gluon MHV amplitude 

We expect that in the light-cone limit (6.3) the five-point correlator (6.1) is related to the five- 
gluon planar MHV amplitude. The latter has been computed to two loops in [31] and it has the 
form 

A 5 = Af [1 + aM^ + a 2 M( 2 ) + 0(a 3 )] , (6.9) 

where A 5 is the tree five-gluon MHV amplitude. Here the one- and two-loop corrections are 
given by the following expressions 

cyclic 

M(2) = \ E [*W< 2) + sj 2 s 15 li 2) + s 12 s 34 s 45 lP] + M 2 (odd) , (6.10) 

cyclic 

where Sij = (pi + Pj) 2 is the invariant mass of the gluons with labels i and j. All gluons are 
considered outgoing, so that their momenta satisfy the relations Y^iPi = an d P 2 = 0. 

Further, I a is a one- loop and I a , l{ , , Ic are two-loop planar scalar integrals depending 
on the gluon momenta, Eqs. (5.12) and (5.13). Switching to dual coordinates, Pi = Xi — Xj+i, 
these integrals can be expressed in terms of the basis integrals (6.8): 

1^=9(1, 2, 3, 4), I® =h(l,3, 2; 1,3, 4), 

I {2) = h(l, 3, 2; 1, 3, 5) , IP = p(2; 1, 3; 4, 5) . (6.11) 

Finally, M odd and M odd in (6.10) stand for the parity-odd contributions to the five-gluon 
MHV amplitude, proportional to the pseudo-scalar e(pi,P2,P3,P4)- The one-loop parity-odd 
contribution M odd is of order 0(e), whereas M odd has a simple pole 1/e. The residue at this pole 
is proportional to the product of the one-loop parity-even and parity-odd parts. As a consequence, 
the parity-odd contribution can be significantly simplified by considering the logarithm of the 
ratio A$/A b , as shown in [31]. Namely, we introduce the following ratio function 

F A = In (A 5 /4 0) ) 2 = aF^ + a 2 F {2) + 0(a 3 e ) , (6.12) 

FJ 1} = 2M« , F {2) = 2M<2) - (M«) 2 . 
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The main advantage of (6.12) is that, unlike the amplitude itself, the parity-odd contributions 
to F A and F A vanish as e — > 0. Replacing M ^ and M^ in (6.12) by their explicit expressions 
(6.10) and going to dual coordinates Pi = Xi — Xj+i, we obtain (up to terms vanishing as e — > 0) 

F ( a ] = 4^242(1,2,3,4) + (cyclic) (6.13) 

Ff = 4^24^(1, 2, 3; 1, 3, 4) + x 4 24 x 2 13 h(l, 2, 4; 2, 3, 4) + 4 4 4 5 4 5 p(l; 2, 5; 3, 4) 

- M 3 xi 4 g(l : 2, 3, 4)] 2 - \x\Ag{l, 2, 3, 4)4 4 4 5 <?(2, 3, 4, 5) 

- |4 3 44^(1, 2, 3, 4)4 5 4 4 <?(3, 4, 5, 1) + (cyclic) . (6.14) 

Here g, h and p stand for the dimensionally regularized integrals defined in (6.8). 

Although the explicit expressions for individual two-loop integrals in (6.14) are quite compli- 

(2) 

cated, their sum F A can be found in a closed form thanks to the remarkable duality between 
MHV amplitudes and light-like Wilson loops [2, 3, 4]. This duality allows us to formulate the 
dual conformal Ward identity for Fa to all loops [6] . For five points it has a unique solution (up 
to an additive constant) which coincides with the BDS ansatz [15]. 

6.2 Duality correlator/amplitude 

Assuming the duality relation (3.19), the correlator on the light cone should be related to the 
scattering amplitude 

lim In (G 5 /GP) = In U,/Af) 2 + 0(e) . (6.15) 

In terms of the functions Fq and Fa introduced in (6.5) and (6.12), this relation implies that up 
to terms vanishing for e — > 

F<P=Fg>, F A *> = F<f, ... (6.16) 

The first of these relations follows immediately from the explicit one-loop expressions (6.13) and 
(6.6). At two loops, F A and F G are given by two seemingly different expressions, Eqs. (6.14) 
and (6.7), respectively. Then, for the duality relation (6.16) to hold, the two-loop integrals in 
(6.14) and (6.7) should satisfy a very non-trivial identity. 

To begin with, let us examine the topology of the Feynman integrals in the two-loop expression 

(1) (1) 

for F A and Fq. By construction, the planar MHV amplitude A 5 receives contributions from 
planar Feynman diagrams only. However, taking the log of the two-loop amplitude (6.9), we find 
that F A contains an admixture of products of one-loop integrals in (6.14). For the correlator the 

(2) 

situation is different. The diagrammatic representation of Fq is shown in Fig. 5. We observe 

(21 

that the expression for F^ involves two-loop integrals of the same type as those contributing to 

(2) 

F A , Eq. (6.14). Closer examination shows, however, that the external points in these integrals, 
xi, . . . ,x$, are not (anti) clock-wise ordered and, therefore, these integrals do not contribute to 

(2) 

the planar amplitude F A . Still, they contribute to the two-loop planar correlator. This is due 
to the different realizations of the planarity condition for correlators and scattering amplitudes. 
For the latter, the planar diagrams have the topology of a disk, while for the former they have 
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/i(l, 2, 3; 1,3, 4) 



/i(l, 2, 4; 2, 3, 4) 



/i(l, 2, 4; 1,3, 4) 
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MM, 5; 1,3,4) 



Ml, 2, 3; 1,4,5) 



Ml, 2, 4; 1,3, 5) 






p(l;2,4;3,5) 



p(l;3,5;2,4) 



p(l;2,5;3,4) 

Figure 5: Two- loop pseudo-conformal integrals of different topologies contributing to the correlator G§, 
Eq. (6.7), and to the amplitude A$, Eq. (6.14). The diagrams with solid lines depict Feynman integrals 
in x— space. The diagrams with dashed lines represent the same integral in the (dual) momentum space 
Pi = %i — %i+i- In the latter case, (— fc) stands for the particle with momentum (—pk)- Thin solid lines 
denote numerators in the x— integral. In momentum space, this numerator is given by the squared sum 
of the momenta flowing through the arrowed dashed lines. 



the topology of a sphere (see Fig. 7 in Appendix A. 3. 2 for an example). Despite of this, the 

(2) (2) 

expected duality relation F y A ' — Fq suggests that the two sums of two-loop integrals of different 
topologies, Eqs. (6.7) and (6.14), are equal to each other. 

We can apply the relation Pi = Xi — Xj + i to rewrite the Feynman integrals contributing to 
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(2) 

Fq as conventional momentum integrals, as shown in Fig. 5. The absence of cyclic ordering of 
the external points Xi leads to an unusual feature of the resulting momentum integrals: They 
involve seven external legs including one pair of legs with opposite momenta, pi and (—pi). 

6.2.1 Soft limit 

There exists a simple way to test the duality relations (6.15) and (6.16) by examining them in the 
special limit x 4 — > x 5 . For the amplitude, this corresponds, e.g., to the soft limit p 4 = x 45 — > 
when one of the external gluon momenta vanishes. In this limit, the dimensionally regularized 
five-gluon MHV amplitude A 5 reduces to the four-gluon amplitude A 4 in such a way that the 
ratio function (6.12) takes the form 

lim F A = In (AJA^Y . (6.17) 

For the correlator (6.1), the limit x 4 — > x^ corresponds to its short distance asymptotic behavior. 
In this limit, we apply the operator product expansion 0(4)0(5) ~ (0 + (4)0 + (5))(9(4) + . . ., 
where (00) is a free scalar field propagator in D = 4 — 2e dimensions, and the dots denote 
subleading terms. In this way, we find from (6.5) that the ratio of the 5-point correlators reduces 
to the ratio of the 4-point ones G(l,2,3,4) = (0(1)0(2)0 (3) 6(A)), 

u ™ Fg = 2 lim ln (?wrA%) ' ( 6 - 18 ) 

where the light-cone limit in the right-hand side corresponds to xf 2 , x 23 , x 2 4 ,x 4l — > 0. 

Let us verify (6.17) using the expression for the two-loop ratio function (6.12). We recall 
that, by definition, the adjacent points in dual space are light-like separated, xf i+1 = 0. In 
the soft limit x 4 — > x^ we find the additional relations x 2 5 — > and x 4l — > 0. Notice that the 
dimensionally regularized g— , h— and p— integrals, Eq. (6.8), remain finite as x 4 — > x^. The 
expression for Fa simplifies because some of these integrals are multiplied by x\ b and x\ x , so they 
do not contribute in the soft limit. Thus, the pentabox integral p drops out: 

lim Fj 2) = 2 xj 3 xl 4 h(l, 3, 2; 1, 3, 4) + 2 x^x 2 13 h(2, 4, 1; 2, 4, 3) - [x 2 13 x 2 24 g(h 2, 3, 4)] 2 . (6.19) 

X4,— >X5 

Using Eq. (5.15), we verify that the expression in the right-hand side of (6.19) coincides with the 
two-loop correction to the ratio of four-gluon amplitudes 2\n(A 4 /A i '), in agreement with (6.17). 
In a similar manner, we examine the ratio function Fq and simplify the expression (6.7) 
in the limit x 4 — > x§. We find that, unlike the case of F A , the pentabox integrals have a 
non-vanishing contribution which can be expressed in terms of h— integrals. For instance, 

~4a 4A P(l; 3, 5; 2, 4) + (cyclic) ^ -6(1, 2, 3, 4) - xi 3 x 2 24 h(l, 3, 2; 1, 3, 4) , (6.20) 

where the notation was introduced for the 'bad' integral 6(1,2,3,4) = x| 4 xf 3 /i(2, 3, 1; 2,4,4), 
which does not have an interpretation in terms of two-loop four-particle planar momentum in- 
tegrals. The same 'bad' integral comes from h(l, 2, 5; 1, 3,4) and p(l;2,4;3, 5) and their cyclic 
images. We find that 6(1,2,3,4) cancels in the sum of all terms in the right-hand side of (6.7) 
as x 4 — > x$, leading to 

lim F G 2) = 2 x\ 3 x 2 2A h(l, 3, 2; 1, 3, 4) + 2 x 4 24 x 2 13 h(2, 4, 1; 2, 4, 3) - [x 2 3 ^ 4 ^(l, 2, 3, 4)] 2 . (6.21) 

X4— >x$ 
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Again, we verify that the expression in the right-hand side coincides with the ratio of four-point 
correlators (6.18). 

Comparing (6.19) and (6.21), we conclude that the two-loop corrections to Fa and F G coincide 
in the soft limit x 4 — >■ x 5 , 

<*) 77,(2) 



lim F ( A > - FX> = . (6.22) 

(2) (2) 

In the next subsection we argue that for arbitrary Xi the difference F A — F G ' is a constant. 
Together with (6.22) this immediately implies that F A = Fq, as announced in (6.16). 

6.2.2 Conformal symmetry and proof of the integral identity 

(2) (2) 

As was already mentioned, the duality relation F A = F G implies an identity between two-loop 
integrals of various topologies. Evaluating the difference F A — F G with the help of (6.7) and 
(6.14) and equating it to zero, we formulate the identity that we expect to find as follows: 

= x 2 13 x 2 24 [x 2 13 h{l, 2, 3; 1, 3, 4) + x 2 24 h(l, 2, 4; 2, 3, 4) - x 2 4 h(l, 2, 4; 1, 3, 4 
+ x\ 3 x\ A x 2 25 [2 h(l, 2, 5; 1, 3, 4) - h(l, 2, 3; 1, 4, 5) - h(l, 2, 4; 1, 3, 5)] 
+ AA [4;p(l;2,5;3,4) - ^ 4 p(l;2,4;3,5) - x 2 35 p(l; 3,5; 2,4 



.?:, 



2 13 x 2 2A g{l, 2, 3, 4)] 2 + \x 2 13 x 2 2i g{l, 2, 3, 4)x 2 24 x 2 3b g(2, 3, 4, 5) 



1~2 



— t;X 



2 ^ l3 x 2 2A g{\, 2, 3, A)x 2 3b x\ i g{l, 3, 4, 5) + (cyclic) . (6.23) 

In the previous subsection we showed that this relation holds in the soft limit. 

The proof of (6.23) goes as follows. The Feynman integrals in (6.23) are regularized dimen- 
sionally with D = 4 — 2e. We will first show that the sum of Feynman integrals in the right-hand 
side of (6.23) is finite as e — > 0. This will allow us to remove the regulator, i.e. to restore the 
D = 4 integration measure. Since all integrals we are dealing with are pseudo-conformal, the 
expression in the right-hand side of (6.23) will thus become an exactly conformally invariant 
function of x±, ... ,x$. Since one cannot construct conformal cross-ratios from the five external 
points with light-like separated neighbors x 2 i+1 = 0, conformal invariance will imply that the 
right-hand side of (6.23) is a constant. We have already seen that the right-hand side of (6.23) 
vanishes in the soft limit, so the constant must be zero. 

The direct way to prove finiteness is to evaluate the pole part of (6.23) by the Mellin-Barnes 
(MB) method (see, e.g., [38]) using the package [39] for the e-expansion and evaluation of the 
representations. At 0(l/e 4 ) through 0(l/e 2 ) the program finds only one-parameter integrals 
which we have analytically evaluated. At 1/e 4 and 1/e 3 one has to show the cancelation of 
rational numbers and simple logarithms, respectively, which was in either case immediate. At 
1/e 2 we needed Landen's identity on dilogarithms, but once again the proof was ultimately 
simple. With our MB representations, the highest integrals contributing to the simple pole in e 
were of dimension four (from the pentabox). The precision of the numerical evaluation used by 
the package [39] was therefore good: For all kinematic points in our (rather general) sample set 
we obtained values like 0.00(2)/e. Within the given accuracy the sum of integrals is seen to be 
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finite. 9 As we have just explained, the finiteness of the integral implies that it equals zero since 
a no n- vanishing constant part is ruled out by the soft limit. 

Alternatively, we can show finiteness analytically by combining the sum of integrals in (6.23) 
into a single integral of the form 

r.h.s. of Eq.(6.23) = / — 2 -^-— 2 , (6.24) 

J V x 10 • • • X 50/V x 10' • • • ^SO'/^OO' 

where, by construction, the polynomial P(xq, xo>,Xi) is invariant under cyclic shifts of the external 
points Xi (with i = 1, . . . , 5) and under the exchange of the integration points xq -h- xqi. Each 
integral in (6.23) gives a contribution to P(xq, xo',Xi) in the form of a product of seven distances 
x\a with various choices of indices % and j. For instance, the third h— integral in the first line of 
(6.23) produces the contribution (— £i 3 £i4^24 x 30 x 50 x 20' x 50')' while the contribution of the cross- 
product of one-loop integrals in the last line of (6.23) looks as (— l^ia^l^u^is^io^lo'^oo')- To 
save space, we do not present the explicit expression for P(xo,xo>;xi). Due to the pseudo- 
conformal property of the integrals in (6.23), the polynomial P(xo,xo',Xi) is covariant under 
conformal transformations, with weight 4 at points x and Xy. Given the weight (—8) of the 
denominator, and if the integration measures can be made four-dimensional, the integral will be 
conformally covariant, as stated above. 

Let us identify the potential divergences of the integral (6.24). They could only come from 
the part of the phase space of xq and xq', in which some of the propagators in the denominator 
of (6.24) vanish simultaneously. Notice that upon the change of variables, p\ = xu, ■ ■ ■ ,P5 — £51 
(with p 2 = 0) supplemented with k = x 10 and k' = x w r, the integral (6.24) can be rewritten as a 
conventional two-loop momentum-space Feynman integral with k and k' being the loop momenta 
and pi,...,p5 defining the external leg momenta. It is well known that infrared divergences 
originate from integration over the loop momenta collinear to the light-like momenta of the 
external legs. For instance, suppose that k^ is collinear either to pi and to p^- To analyze 
the divergences, it is convenient to employ the so-called Sudakov decomposition of the loop 
momentum, 



x 
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apt + prf + k^_ , (6.25) 



where a, (5 are scalar variables and k^_ are two-dimensional transverse momenta, (k±pi) 
{k±P2) = 0. In terms of Sudakov's variables, the xq— integral in (6.24) takes the form 



r d 4 x r x 2 13 dad/3d 2 k ± [...] 

J xMoxlo "- 1 ~ J (rt&P - kl)(xUl - a)(5 + k\){x\ 



(xj 3 a(3 - k 2 ± )(xl 3 (l - a)(3 + kl)(xj 3 a(l - 0) + k\) ' 



(6.26) 



where [. . .] denotes the remaining terms in the right-hand side of (6.24). 

It is easy to see that the integral (6.26) develops logarithmic divergences originating from the 
integration region 

x 2 13 a/3 -k\ = 0(p 2 ) , k± = O(p) (6.27) 

with p — > 0. Depending on the hierarchy between a and /3 we can distinguish three subregions: 
For a,(3 = O(p) we have k^ = O(p), for a = 0(p°), (3 = 0(p 2 ) we have k^ = apl + O(p) and, 



9 An analytical proof of the absence of the simple pole, too, should be feasible given the relatively low dimen- 
sionality of the MB-integrals. 
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finally, for a = 0(p 2 ), /3 = 0(p°) we have k^ = ftpi?, + 0(p). In terms of the dual x— variables, 
this corresponds to the limit where the integration point Xq approaches either the external point, 
Xq — ¥ xi, or one of the light-like segments [xi,x 2 ] and [#2, £3], respectively. Due to the cyclic 
symmetry of the integral (6.24), divergences are also produced when xq approaches the other 
cusp points Xi and the light-like segments [x^Xj+i]. The same analysis applies to the integral 
with respect to xq>. 

Then, for the integral (6.26) to be finite, the expression inside [. . .] should be finite and, 
in addition, it should vanish in the region (6.27). This leads in its turn to the condition for 
the polynomial P(xa,xy,Xj) in the right-hand side of (6.24) to vanish sufficiently fast when xq 
and/or xq> approach one the potentially dangerous regions explained above. We return to (6.24) 
and verify that the polynomial P(xo,xq'] xi) indeed satisfies this condition. Thus, the integral 
(6.24) remains finite as e — > 0. Together with conformal invariance this immediately leads to the 
identity (6.23). 

7 Conclusions 

The main result of this paper is the observation that computing the loop corrections to the corre- 
lators of protected operators by the Lagrangian insertion method, and taking the light-cone limit, 
gives us the integrands of the MHV gluon scattering amplitudes in the dual momentum space. 
These integrands can be evaluated directly in four dimensions since they are finite and explicitly 
(dual) conformally covariant. The divergences appear in the integrals over the insertions points 
for the correlator (UV divergences) or, equivalently, over the loop momenta for the amplitude (IR 
divergences). In order to achieve an exact matching of the two objects, we need to dimensionally 
regularize the theory, not as one would do for a correlator computation in coordinate space, but 
in the way that is natural for the dual theory in momentum space. We call this dual infrared 
regularization. 

At present we have no understanding of the field-theory mechanism responsible for this sur- 
prising relation between two seemingly very different objects. In fact, the surprise is not total, 
since we already know that MHV amplitudes are dual to light-like Wilson loops. The latter, as 
shown in [28], are intimately related to correlators in the light-cone limit. 

Further tests are needed to confirm our conjecture. In particular, a crucial test will be the 
calculation of the six-point correlator at two loops and its comparison to the six-gluon amplitude. 
This is similar to the test of the Wilson loop/amplitude duality at six points and two loops [7, 8] . 
Before it, one could suspect that the matching of the four- and five-point objects was simply due 
to (dual) conformal symmetry, which fixes their form completely. As we saw in our five-point 
two-loop correlator computation in Sect. 6, conformal symmetry plays an important role here 
too. However, at six points it will not be sufficient to explain the matching, if it is confirmed. 

The Lagrangian insertion procedure described here can give us a simple way to compute the 
integrands that determine the higher-loop corrections to amplitudes. We would like to emphasize 
that our one- and two-loop correlator computations involve a handful of Feynman graphs. This 
is in striking contrast to the huge number of graphs needed, if one would attempt to do a 
straightforward amplitude calculation (not using unitarity methods). 

The fact that the integrand of the amplitude can be obtained from a tree-level computation 
raises the hope that one might be able to profit from some hidden enhanced symmetry of the 
tree-level correlator. This is motivated by the analogy with the tree-level scattering superampli- 
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tudes, which are known to have a dynamical symmetry called dual conformal symmetry [40]. If 
something similar exists also for the correlators, it might eventually lead us to the integrability 
of the loop amplitudes. 

A related issue is the observation that the correlators we have discussed enjoy the full super- 
conformal symmetry of the M = 4 SYM theory. In this paper we have explored the lowest scalar 
components of the supersymmetric correlators, and we have shown that they match the gluon 
MHV amplitudes. The question arises if the correlator defined in superspace has some relation 
to the superamplitude, and thus possibly to non-MHV amplitudes. 
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Appendices 

A Lagrangian insertion procedure in J\f = 2 harmonic su- 
perspace 

In this appendix we give a very brief overview of harmonic superspace, in particular of the 
Feynman rules we need here. More details can be found in [34, 35, 41]. We then summarize 
the Lagrangian insertion procedure, developed in [42, 23] for the case of four-point correlators. 
Finally, we explain how to adapt this procedure to five (and more) points. 

A.l M = 4 SYM in M = 2 terms 

The basic ingredients of the M = 4 SYM theory in M = 2 terms are the hypermultiplet and the 
super- Yang-Mills (SYM, or vector) multiplet. The M = 4 SYM action consists of two terms: 

Sn=4 SYM = Sn=2 SYM + ^HM . (A.l) 

Below we give a short description of each multiplet and its action. 

A. 1.1 M = 2 hypermultiplet in harmonic superspace 

The M = 2 massless matter (or hyper) multiplet consists of an R-symmetry SU(2) doublet of 
complex scalars 4> l {x) (with i — 1,2) and of two Majorana spinors and SU(2) sunglets ip a (x), 
ipa{x) and K a (x), Ra{x). Their supersymmetry transformations close only on shell. Going off 
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shell requires the introduction of an infinite set of auxiliary fields [34]. This is achieved by 
extending the space-time by two extra compact dimensions in the form of a sphere S* 2 . The 
latter is described in terms of harmonic variables u ±l which form a matrix of SU(2), 

|| u || 6 SU(2) : u +i uT = 1 , ~i = u r = eij u- j , ei2 = -e 12 = 1 , (A.2) 

and parametrise the sphere S 2 ~ SU(2)/U(1). A harmonic function /^(m 1 * 1 ) of U(l) charge q 
is a function of u ±l invariant under the action of the group SU(2) (which rotates the index % of 
u ±l ) and homogeneous of degree q under the action of the group U(l) (which rotates the index 
± of u ±l ). Such functions have infinite harmonic expansions on S 2 whose coefficients are SU(2) 
tensors (multispinors). 

In this framework the hypermultiplet is described by a harmonic superfield q + (x,9,9,u) of 
U(l) charge +1 satisfying the Grassmann (or G-)analyticity constraints 

DU + = D+q + = , (A.3) 

where 

D+ = D\uj , D+ = Dlut (A.4) 

and D l a , D\ are the usual supersymmetric spinor derivatives. These constraints can be solved 
explicitly in the G-analytic basis in superspace 

xf = x ac * - Aid a H^utu- , B% & = ute i aA (A.5) 

where x aa = x^a^ and (ij) means weighted symmetrization. In this basis q + becomes a function 
of 6 + ,6 + only, i.e., a G-analytic superfield q + (xA, + , 6 + ,u). 

As mentioned earlier, the M = 2 supermultiplet can exist off shell because an infinite number 
of auxiliary fields (coming from the harmonic expansion on S 2 ) are present. On shell these 
auxiliary fields are eliminated by the harmonic (or H-)analyticity condition (equation of motion) 

D ++ q + = . (A.6) 

Here D ++ is the harmonic derivative on S* 2 (the raising operator of the group SU(2) realized 
on the U(l) charges, D ++ u + = 0, D ++ u~ = u + ). In the G-analytic basis (A.5) it becomes a 
supercovariant operator involving space-time derivatives: 

D ++ = u +l —- - 4i9 +a 9 +A -^- . (A.7) 

It is then easy to show that the free on-shell hypermultiplet becomes an "ultrashort" superfield: 

q + (x A ,9 + ,9 + ,u) = <P\x A )uj + 9 +a ^ a (x A ) + 9tK«(x A )+Ai9 + a»9 + d ll <P\x A )u- , (A.8) 
where the physical scalars 0* and spinors ip a , R a satisfy their massless field equations D0*(x) = 

The equation of motion (A.6) can be derived from an action given by an integral over the 
G-analytic superspace: 



Shm = - 2 / dud 4 x A d 2 9 + al 2 9 + Tr (q + D ++ q + ) . 



(A.9) 



Here q + (x A , 9 + ,9 + , u) is the conjugate of q + . The conjugation ~ combines usual complex conju- 
gation with the antipodal map on S 2 in a way to preserve G-analyticity. This action is real (with 
respect to the ~ conjugation) which can be seen by integrating D ++ by parts. In this sense the 
action (A.9) resembles the Dirac action for fermions, although the superfield q + is bosonic. 
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A. 1.2 J\f = 2 SYM multiplet in harmonic superspace 

The M = 2 SYM gauge potential is introduced by covariantizing the action (A. 9) with respect 
to a Yang-Mills group with G-analytic parameters X(x A ,9 + , 9 + ,u). To this end one replaces the 
harmonic derivative in (A. 9) by the following covariant one: 

D ++ ->• D ++ + igV ++ (x A , 9 + , 9 + , u) , (A.10) 

where g is the gauge coupling constant. The gauge potential is described by a real (V ++ = V ++ ) 
G-analytic superfield of charge +2 (equal to the charge of D ++ ). The matter and gauge superfields 
are subject to the usual gauge transformations: 

q+' = e l9X q + , V ++ ' = --e igX D ++ e~ i9X + e iaX V ++ e- l9X , (A. 11) 

9 

so that the covariantized action (A. 9) 

Shm/sym = -2 J dud A x A d 2 9 + d 2 e + Ti{q + D ++ q + + ig q + V ++ q + ) (A. 12) 

is indeed gauge invariant. 

In the non-supersymmetric Wess-Zumino gauge the gauge potential has the component ex- 
pansion 

V++(x A , 6 + , 6 + , u) = -2i9 + a»6 + A,(x A ) - iV2($ + ) 2 <p(x A ) + iy/2(0+) 2 ip(x A ) (A.13) 

+ 4(6+) 2 e+<*X i a (x A )u- - A{9 + ) 2 9t\ & \x A )u- + 3(6+) 2 (d+) 2 YV(x A )u;uJ , 

containing the fields of the M = 2 off-shell vector multiplet: the gauge field A M , the complex 
physical scalar tp, the doublet of Majorana gluinos X l a , X at and the triplet of real auxiliary fields 

The gauge invariant action for V ++ can be written down either directly in terms of V ++ 
[41, 35], or in terms of the chiral superfield strength W(xL,6 ia ) (or its conjugate antichiral 
W(x R: 9 idt )): 

Sn=2 sym = ^2 / d4xLd " 9 TtW2 = Y^[ faRffi Tr W 2 , (A. 14) 

where 

x ™ = x «« - 2i9 ia 9? (A. 15) 

are the space-time coordinates in the chiral basis, and xr = xl are the antichiral ones. In a topo- 
logically trivial background these two forms are equivalent (up to a total space-time derivative), 
due to the Bianchi identity D ia D{W = D^D&W. 

Unlike the G-analytic potential V ++ , the field strength W(xl,9) is a chiral superfield which 
does not depend on the harmonic variable m ± . It can be expressed as a power series in V ++ , 
involving multiple harmonic integrals [43] 

y -W^g/*»---^ ( gSSg:::SS ) ' ^ 
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where (w + -u + ) = u +l €ijU 1 3 , etc. In terms of component fields we have (in the Abelian case) 
W = <p(x) + e ai X ai (x) + e ai eP j (e i:j (a^UF^(x) + e^Yij) + i{9^d ad X\x) + 9 4 D^(x) . (A.I7) 



With the hypermultiplet matter in the adjoint representation of the gauge group, the two 
actions (A. 12) and (A. 14) describe the M = 4 SYM theory, 

Sn=4 SYM = Sn=2 SYM + Shm/SYM • (A. 18) 

As mentioned earlier, the main advantage of the M = 2 harmonic superspace formulation is 
the possibility to quantize the theory in a straightforward and manifestly M = 2 supersymmetric 
way [41]. Further, compared to the M = 1 chiral matter superfields, the M = 2 hypermultiplet 
composite operators O = Tr(g + ) fc , etc., need no covariantization, hence no presence of the gauge 
superfield in the definition of the correlators (O . . .0). The hypermultiplet matter interacts with 
the gauge sector only through a single trilinear vertex, which considerably simplifies the Feynman 
diagrams. The true non- Abelian nature of the theory is encoded in the gauge self-interactions 
(as well as in the ghost sector, but we do not need it here). 

A. 2 Feynman rules 

In this section we give a subset of the Feynman rules for the combined M = 2 matter+gauge 
system which are sufficient for the one- and two-loop calculations we do in this paper (the 
complete set can be found in [41, 35]). 

The hypermultiplet propagator is the solution to the Green's function equation 

Df + (q + (l)q+(2)} = 8\x A1 - x A2 )5 2 (0+ - 6})5 2 (6+ - 9t)S( Ul , u 2 ) , (A.19) 

and is given by 



(e(l)% + (2)> = 4^r*«* = > (A.20) 



(12) '-" -''■ 

-12 

Here a, b are color indices and 

(12) = -(21) = ut%jUp (A.21) 

is a shorthand for the SU(2) invariant but U(l) covariant contraction of the two harmonics. The 
coordinate difference 

2? - - 

x1 2 = x» A1 - x% + [(1-2)0+^0+ + (2-1)0+0*0+ + 0+o"0 2 + + 2 V0 + ] , (A.22) 

where, e.g., (1 _ 2) = Ui' l tijU 2 \ is invariant under the Poincare supersymmetry transformations 
in the G-analytic basis (A.5): 10 

S Q xf = -4ni"(e iQ 0+ d + 0+V d ) , 5 Q 6 +aA = u+e ia ' d , 5 Q uf = . (A.23) 

Setting the Grassmann variables in (A.20) to zero, we find the propagator for the physical 
scalars (f> l a (x) projected with harmonics (recall (A. 8)): 

(e(lk + (2)>.=o = (<>„(*!) K-^(:r 2 )> = j^Sa, ■ (A.24) 

47T X 12 



10 To check this one makes use of the harmonic cyclic identity (1 2)1 + (21)1 + (11 )2 = and of the defining 
property (II - ) = 1 (sec (A. 2)). 
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The gauge field (gluon) propagator depends on the gauge we have chosen. In our loop 
calculations we will only need the propagator (W(l)V ++ (2)), having one chiral end (the field 
strength W(xli,9i)) and one G-analytic end (the SYM potential V ++ (xa2, &t ■> ®t-> u 2))- It is 
independent of the gauge and has the following form 



1,0 2,6 



(W a (l)vr(2)) = -^gr(M 2 = #^00000 (A.25) 



It involves the coordinate differences 

xt = xz - x a A - ^e? el" , e? 2 = v+e? - e+ a (A.26) 

and is invariant under the Poincare supersymmetry transformations in the chiral basis for xli, 

8 Q xf = -4ie ai 4 , S Q e ia = e ta , (A.27) 

and (A. 23) in the G-analytic basis for xa2- 

Finally, the only interaction vertex we shall need here is the gluon-to-matter coupling which 
can be read off from the covariantized hypermultiplet action (A. 12): 

vr 

= igfbac I d 4 x du d 4 9^ (A.2£ 

Qa 

^— 



The SYM action (A. 14), (A. 16) is non-polynomial in V ++ , so it contains infinitely many vertices 
(but, of course, becomes polynomial in the Wess-Zumino gauge (A. 13)). At the one- and two- 
loop levels that we are studying in this paper only the cubic non-Abelian vertex can appear, but 
later on we shall see that all such graphs are irrelevant. 

A. 3 The insertion procedure: Four points 

Here we illustrate the M = 2 Lagrangian insertion procedure described in Sect. 3, on the example 
of the two-loop four-point correlator G 4 (5.1) of hypermultiplet bilinears O = O d + = §+ =0 with O = 
Tr(q + ) 2 and O = Tr(g + ) 2 . In particular, following [23] we explain the important role of M = 2 
superconformal symmetry, which allows us to drastically simplify the two-loop calculations. 

The two-loop (order g A ) corrections to the four-point correlator G4 can be obtained by a 
double insertion of the M = 2 SYM chiral action 

Sm=2 sym = / d 4 :rd 4 # L(x L , 6), L = l^Tr(W(x L , 9)f . (A.29) 

To see this, we first rescale the gauge potential V ++ — >■ g~ 1 V ++ , with the effect that the coupling 
only appears in front of the SYM action (A. 14), but not inside the field strength W (A. 16). It 
also drops out of the gauge/matter coupling (A. 12), thus (A. 28) loses the explicit g. As another 
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consequence, the gauge propagator (i.e. the inverse of the gauge kinetic term) is scaled up by g 2 , 
which introduces a factor of g 2 in the right-hand side of (A. 25). 
Let us write out the perturbative expansion of 4 : 

G A = of + g 2 G { P + g 4 Of + . . . , (A.30) 

so that 

On the other hand, by considering the effect of the differentiation on the path integral we find 

\ gA {w) ° 4 = ^/^o^o (^(0)0(1)0(2)0(3)0(4)) (A.32) 

-- f d 4 x d 4 9 d 4 x ,d 4 9 , (L(0)L(0')O(l)6(2)O(3)O(4)) . 

The left-hand side in the last formula starts at 0(g 4 ), which leads to a puzzle: The perturbative 
expansion of the first term on the right-hand side starts at 0(g 2 ), while the second term seemingly 
starts at 0(g 4 ). The only way to produce a compensating 0(g 2 ) contribution from there is to 
insert both L(0),L(0') into one gluon line. This means inserting the chiral-to-chiral propagator 
(it contains just the propagator of the auxiliary field Y in (A. 17)) 

(W(O)W(O')) = 2tg 2 5 4 (x L0 - x L0 ,)5 4 (9 - O v ) (A.33) 

into that gluon line. Upon performing the chiral superspace integration over point 0', the 0{g 2 ) 
contribution from the single insertion term is identically canceled. This remains true in general: 
00' contact terms from the double insertion term identically cancel the single insertion term. 
Returning to 0(g 4 ): 

Of = ~ j d 4 x d 4 # dVdX' (L(0)L(0')O(l)O(2)O(3)O(4)>; e 4 g (A.34) 

where the superscript indicates that the 00' contact terms are to be omitted. For the class of 
graphs we find below this simply means not to contract W(0) and W(0'). We stress that this 
(4 + 2)— point correlator is once again at Born level; it comes with g 4 because there are four 
Yang-Mills propagators and two explicit factors g~ 2 from the Lagrangian insertions. Below we 
show that this correlator with insertions has to be nilpotent, i.e. proportional to 9 s . 

A. 3.1 Structure of the nilpotent superconformal covariant 

The most important feature of the new six-point correlator is its superconformal symmetry. 
Indeed, it involves gauge-invariant composite operators, O, O and L. As explained in Sect. 3, 
all of these operators are particular projections of the M = 4 half-BPS protected operator in 
the 20' of SU(4). As such, they need no renormalization and have well-defined superconformal 
properties. As long as we keep the end points in this correlator apart, nothing can break the 
M = 2 superconformal symmetry of the theory. This imposes rather strong constraints on the 
general form of the correlator. 



31 



The M = 2 superconformal algebra has an SU(2) x U(l) automorphism group (R symmetry). 
The U(l) factor (to be distinguished form the harmonic U(l) C SU(2)) acts only on the odd 
superspace variables, R[9] = 1/2, R[9] = -1/2. From the SYM action (A.29) we deduce R[L] = 2 
and i2[W] = 1. At the same time, the hypermultiplets q + and q + have no R charge, as follows 
from the action (A. 9), but they carry harmonic charge 1. This implies that the six-point correlator 
carries harmonic U(l) charges 2 at points 1 to 4, and a total R charge 4 at points and 0'. Since 
the chiral 9 a (at both the insertion and external points) are the only superspace coordinates with 
positive R charge 1/2, we conclude that the correlator can be written in the factorized form 

(LLOOOO) = 0(x, 9, u) x f(x, u) + fl-terms , (A.35) 

where G is a particular nilpotent six-point superconformal covariant, homogeneous in 9 of degree 
8, and thus carrying the whole R charge. The antichiral odd variables 9 can only come from the 
external points, but we are ultimately interested only in the lowest components O = O e + = §+ =0 , so 
we can ignore the 9 terms in (A.35). So, the essential information about the six-point correlator 
(A.35) is contained in the function f(x,u) without R charge. Below we will show that this 
function is in fact harmonic independent, the harmonic U(l) charge being carried by 0. 

The structure of the nilpotent covariant is determined by superconformal symmetry com- 
bined with the G-analytic nature of the four external points and the chiral nature of the two 
insertion points. In addition to Poincare supersymmetry (parameters e l a , ef) we need to consider 
special conformal supersymmetry (parameters rf a , fjf). In the chiral basis (A. 15) we have 

Sxf = -Ai9 ia ef - Ai9 ia x1%i , 59 ia = e ia + xf ffy + 0{9 2 ) , (A.36) 

while in the G-analytic basis (A. 5) we find 



?',; 



Sxf = -4iu-{e ia 9+« + 9 +a ^) + 4i{xf9 + ^ - x^9 + %)^ 

S9 +a = ufe ia + x a /r ] lui + 0(9 2 ), 59 + « = u+e* - xf^uf + 0(9 2 ) , 
5uf = Ai(9 +a rii + f]i9 +& )u+u- , 8u~ = . (A.37) 

Now, the covariant is homogeneous in 9 of degree 8. The inhomogeneous part of 69, 

6d ia = e ia + x ^-i^ ^ 50 +a = u + £ to + x fffy+ , (A.38) 

would lower this degree, unless we find combinations of 0's which are invariant (to lowest order 
in 9, 9) under M = 2 superconformal supersymmetry. Such combinations are 

£ra = Pra ~ &ra , T = 1, . . . , 4 (A. 39) 

with 

Pra = (9? ~ 9 i Q U+ i ) a ( X r0)*aXrt , °ra = (9+ - 0J,U^) Q '(Xr<y) a&X$ . (A.40) 

Their total number is 8, and we wish to construct the nilpotent covariant of degree 8. We 
conclude that the leading term of must involve all of the variables £ r : 

Q=mm- (A.4i) 

The aim of our two-loop calculation is to determine the factor f(x,u) in the six-point corre- 
lator (A.35). Since we are only interested in the lowest component of the four-point correlator 
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(OOOO) = (OOOO) e + = g+ =0 , we can set all the external 9s to zero, 6+ = 0, r = 1, . . . , 4. In this 
case is rather simple [23]: 

ei* + =o = eX 3 0,) T 2 , (A.42) 

1 ir=l X rO X rO' 

where 

U = (12) 2 (34) 2 x 2 4 x 2 3 + (U) 2 (23) 2 x 2 12 x 2 u + (12) (23) (34) (41) [x 2 13 x 2 24 - x 2 12 x 2 u - x 2 u x 2 23 ] . (A.43) 

Finally, substituting everything into the double-insertion formula (A. 32) and performing the 
trivial chiral integrations over #o,o', we obtain the two-loop correlator 

(OOOO) = (OOOd) e+=s+=0 = K f d ^odV (x 2 q/)2 /(a , u) {AM) 

J llr=l ^rO^rO' 

Notice the characteristic presence of the polynomial prefactor 1Z. As shown in [37, 44, 45] , this 
factorization of the loop corrections is a universal feature, called "partial non-renormalization" . 

A. 3. 2 Feynman graphs. Harmonic analyticity 

Now, the practical question is how to compute f(x,u) from the corresponding set of two-loop 
Feynman diagrams. It turns out that instead of setting 6 + = 0, as required in the final expression 
(A. 44), it is much more convenient to do the computations with # ,o' — 0. The knowledge of the 
complete (A. 41) allows us to easily switch from one of these forms to the other. The new form 
of is even simpler, yielding 

(LLOOOd) do ol=§+=0 = (x 2 w r fl -W- /(*> «) • (A.45) 

r=1 X rO X rO' 

Then it is clear that in working out the expressions for the various Feynman graphs we can 
concentrate only on the terms with the maximal number of external 9 + . In particular, at order 
g 4 this choice removes all graphs which contain non-Abelian interaction vertices. For example, the 
non-Abelian gluon subgraph in Fig. 6(c) vanishes because it has two chiral ends at the insertion 
points and 0' and one G- analytic end (the gluon without insertion); after setting o ,o' = we 
are left with too few chiral 6 + s at the G-analytic gluon end to supply the required R charge 2. 
Similarly, the block in Fig. 6(d) has three chiral ends (in fact, only two, points 0' and 0" should be 
identified) and two G-analytic ends; once again, the G-analytic 9 + s cannot provide the required 
R charge 3. The same applies to the block in Fig. 6(e). 

As a result of all these simplifications our task is reduced to listing all tree level Feynman 
graphs made out of the two building blocks T and TT in Figs. 6(a) and (b). They are calculated 
with the Feynman rules from Sect. A. 2. Although these blocks contain gauge/matter interaction 
vertices (A. 28) and hence space-time integrals, the latter are easily done using the identities 

°. I "/^V = 4^ • 3M I ^V = -<**'-&&- , (A.46) 

J x 13" f '23 J '03 x 12 x 10 J x 13 x 23 x 03 X 12 J '10 X 20 

producing very simple rational space-time functions [42, 23]: 

<&(l)W 6 (0) g +(2)> = -p^- [(21-)p 2 + (12-)^ - 2( PlP2 )j (A.47) 

{ZTT) X 12 I J 
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Figure 6: Building blocks of the Feynman graphs. 
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(A.48) 
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where p and a were defined in (A. 40). 

Notice the characteristic presence of negative- charged harmonics in both expressions (A. 47) 
and (A.48). This has to do with the important issue of harmonic analyticity [46]. In an inter- 
acting theory the hypermultiplet satisfies its equation of motion (A. 6) with a covariant harmonic 
derivative, D ++ q + + ig[V ++ ,q + ] = 0. The gauge-invariant composite operators O k = Tr(g + ) fc 
satisfy the same equation with a flat harmonic derivative, D ++ Ok = 0. As explained above, the 
harmonic derivative is the raising operator of £77(2). So, O k corresponds to the highest- weight 
state of an SU(2) irrep of weight k (a (k + l)-plet). 

In practice, this means that the n— point correlator is annihilated by the harmonic derivative 
D ++ at each point, 

D+ + (O k ...O k ) = 0, r = l,...,n. (A.49) 

Since D ++ u + = and D ++ u~ = u + , this implies that the correlator is a polynomial in w+, 
r = 1, . . . , n, homogeneous of degree k at each point, and no dependence on u~ is allowed. 

Clearly, the expressions for the building blocks T (A. 47) and TT (A.48) are not harmonic 
analytic because of the presence of u~^ and u^. This, however, is not a problem: The various 
building blocks or even complete Feynman graphs are not expected to be harmonic analytic, 
much like they are not conformal and gauge invariants. It is only the sum of all graphs that has 
these properties. Indeed, it can be shown that by summing up all graphs made from the T and 
TT blocks, all negative- charged harmonics drop out. To see this one uses the harmonic cyclic 
identity, e.g., 



(12)(1~2- 



;i2-)(l-2) 



uf l Ui j ut k u 



"1 "1 



(dkCji - e u e jk ) = (11 )(22~ 



(A.50) 



as a consequence of the defining property u 



+i 



u. 



1 (A. 2). In practice, the use of the cyclic 



identity is cumbersome when there are w's from too many different points. But we can do better, 
by completely sidestepping this issue. 

We can profit from the expected harmonic analyticity of the final result to greatly simplify 
our graph calculations. Let us come back to the correlator of four operators = 02 with two 
Lagrangian insertions, calculated at 9q = 9q> = 0, see (A. 45). The Lagrangian has no harmonic 
U(l) charge, hence (LLOOOO) should have charges +2 at each external point. From (A. 45) 
we see that the nilpotent factor already carries the necessary charges, thus making the function 
f(x, u) chargeless. Harmonic analyticity then implies that this function is harmonic independent. 
This allows us to compute the correlator (A. 45) with all four harmonic variables identified, 



l/,7 



U% 



Uo 



u,-. 



(A.51) 
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This simple trick eliminates a number of irrelevant Feynman graphs, namely, all graphs with 
at least one free hypermultiplet line, since the hypermultiplet propagator (A. 24) vanishes if 
uf = v,2- Among them we find the graphs with the blocks from Fig. 6(c)-(e), for which we 
already gave a different reason why they do not contribute. In addition, the identification of 
harmonics eliminates the graphs with TT blocks. This leaves only the three graphs shown in 
Fig. 7 (plus point permutations): 






Figure 7: Graphs surviving the identification of harmonics. 

Notice that despite the appearance, the third graph is planar. This has to do with the fact that 
unlike the Green's functions of elementary fields, which are planar on a disk, the correlators 
(Green's functions) of composite operators have the topology of a sphere. This peculiar property 
of the correlators was mentioned in Sect. 6.2. 

Each of the graphs in Fig. 7 is a product of T blocks (A. 47) evaluated at 8q = 9qi = 0. Thus, 
the calculation of the harmonic-independent function f(x) in (A. 45) is reduced to elementary 
algebra. Then f(x) is substituted in (A. 44), leading to the final result (5.4) from [23]. 

This procedure is a very efficient tool for calculating not only four-point, but also n-point 
correlators of half-BPS operators made of hypermultiplets. In the next subsection we show how 
to adapt the procedure to five points. 

A. 4 The insertion procedure: Five points 

Here we evaluate the correlator 

G 5+2 = (1,(0) L(O') 0(1) 0(2) 0(3) 0(4) 6(5)) , (A.52) 

with O = Tr(g+) 2 ,6 = Tr(g+) 2 ,6 = 2Tr(g+g+),L = l/(2g 2 )Tr(W) 2 , in order to study the 
two-loop corrections to the correlator 

G 5 = (0(1) 6(2) 0(3) 6(4) 6(5)) , (A.53) 

discussed in Sect. 6. We will use the technology developed in [23] and reviewed in Appendix A. 3 
for the four-point case G 4 = (0(1) 6(2) 0(3) 6(4)). 

As before, we compute G§ + 2 in D = 4. It remains finite in the light-cone limit 



x 



12 



X 



23 



X 



:-!! 



X 



15 



X 



51 



->0. 



(A.54) 



The divergences of G5 in this limit arise from the integration over the insertion points Xq, Xq> 
which we will eventually regularize by the IR prescription of changing only the dimension of the 
integration measure at these points. 
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Repeating the superconformal argument of Appendix A. 3.1, we can claim that the leading 
term in the 9 expansion of the correlator £5+2 will factorize into a nilpotent covariant ~ 9 8 
and some function f(x,u) of the bosonic variables, see (A. 35). Once again, O will depend only 
on the invariant variables £ r (A. 39) (with r = 1, ... ,5). This time, however, we have 10 such 
variables, while the degree of homogeneity of is still 8. Unlike the four-point case, where the 
nilpotent structure (A. 41) was unique, now we can have two different choices: 



=5/3 . 



^45 = ^^f&d^/^M, (A.55) 



and their point permutations. A graph calculation is needed to fix the coefficient functions 
f(x,u). The knowledge that only A r , B rs can occur is very useful, though, if combined with the 
property of harmonic analyticity, namely the fact that the full gauge invariant correlator only 
depends on u* and that it does so in a polynomial fashion (see Appendix A. 3. 2). 

Let us consider the covariants A r , B rs at 9 = 9 i = 0, in other words, only the terms involving 
9+ from the outer points. Now, 0+ carries U(l) charge 1 at point r, and so does u+ . Hence the 
spinor part of the covariant A^ has charge 2 at points 1,2,3,4, whereas it is chargeless at point 5. 
Since all five outer operators O, O, O carry charge 2, the coefficient function f${x, u) has to have 
charge 2 at point 5 and zero at all other points. To be harmonic analytic it must be an SU(2) 
invariant polynomial in the w+ with the correct charges. The only such invariant is trivially zero, 
(55) = e^u^u^ = 0. This argument rules out all A r covariants. 

Further, in B45 the odd variables £ carry charge 2 at points 1,2,3 and charge 1 at points 4,5. 
Due to harmonic analyticity, the harmonic dependence of the bosonic factor can only be of the 
form f^(x,u) = (45)/f 5 (x). Thus, the sum of the contributions of the graphs to this covariant 
will be 

5 45 = £ il e 3 U e 5/3 (45) fti(x) . (A.56) 

Individual graphs do contain non-analytic terms, but we need not go through all the details 
of how they cancel out. Instead, we can apply the powerful trick of identifying the harmonic 
variables, as we did in the four-point case in Appendix A. 3. 2. 

For each covariant B rs (with r,s = 1, . . . , 5) we know that the harmonic dependence of the 
coefficient function will eventually be given by just (rs). This result clearly does not change if 
all harmonics, except for uf, are put equal to uf. As a convention, if r < s we will keep u^ 
aside and identify all other harmonics with uf { . We will obtain the correct result if this is done 
consistently for any contribution to the spinor structure pertaining to the given covariant (A.56). 

This manoeuvre drastically simplifies the use of the harmonic cyclic identity. For example 
equation (A. 50) reduces to — (—1) = 1 if both uf i: uf { are sent to uf { . But there are more 
far-reaching consequences for the supergraphs at 9q = 9q> = 0. 

• Any diagram with more than one free line (i.e. a hypermultiplet propagator between two 
outer points) is put to zero. Spinor structures relating to the A r covariants may be discarded 
immediately. For the B rs type contributions the suggested identification of the harmonics 
will send at least one of the numerators of the free lines to (ss) = 0. 

• As a consequence, we only need to take into account graphs built out of four T blocks and 
one free hypermultiplet line. 
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• The light-cone limit singles out diagrams in which the gluon lines connect to the pentagon 
frame 123451 of matter lines, because the T blocks with outer ends r, s have an explicit 
propagator factor l/xf s . This eliminates graphs with disconnected matter frames, like the 
third graph in Fig. 7, since they lack the required light-cone singularity. For the same 
reason, graphs with a connected "zigzag" frame like in Fig. 1(b) are not allowed. 

All of these simplifications leave us with a very small number of graphs shown in Fig. 8 (notice 
that the middle graph is planar, like the two others). 






Figure 8: Graphs contributing to 6*5+2. The permutations of the insertion points 0, 0' and the cyclic 
permutations of the external points should be added. 

The resulting algorithm for evaluating these graphs is as follows: 

• For any of the diagrams in this set multiply out the terms of the four T blocks and classify 
them according to the various 8^ structures. Discard terms with four (6^) 2 (A type) but 
keep those related to the B type covariants. Sum over all graphs. 

• In each term identify the harmonics according to the spinor structure; by way of example 



for 



?)W)W) 2 0£0J? we put u± = 4 = 4 = 4. 

may vanish in doing so, in other cases the harmonics will reduce to the simple factor (45). 



The contributions of some diagrams 



• Reconstruct the full covariants. The result of the procedure is the entire leading term of the 
correlator G$ + 2 in which the harmonics are not identified any longer. We can now switch 
to the "opposite end" of the expression by putting all 0+ = 0, while restoring the spinors 
from the insertion points. We observe that this step produces (12)(23)(34)(45)(51) 9q0q, g , 
i.e. in the light-cone limit no other SU(2) channel is present. 

The reconstruction is in fact elementary. From the definitions (A. 39), (A. 40) we find 
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(A.57) 



By this formula we can unambiguously upgrade every 6+ from the graph calculation to the 
invariant combination £ r . Next we note that 
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where 

X£ = ^6?<, te{0,0'}. (A.59) 



x Q0 , 



In order to complete the task we must collect the w-projected A variables into 8q Oq, and a harmonic 
factor. To this end we use the identity 



\1+ \2+ x3+ \4+ 
Od 0/3 O7 05 Q2 



wKt Kt 3 ) A o- -^K)(^Kt-wKt) ( A - 6 °) 



^2 (^^(14)(23) -e^e^. (13)(24)) 9%, 



4 ( x oo - 

and its special cases where some points coincide. Beyond the factor (rs) in the coefficient function 
of the B rs covariants, the conversion to 0^8$, produces four further harmonic factors (ij) which 
now carry the remaining U(l) charges. The harmonic dependence remains manifestly analytic. 
Out of the w+ from the five outer points one can construct six independent polynomials carrying 
charge 2 at every point. We stress that it is not obvious that the sum of graphs produces only 
one channel in the light-cone limit, namely (12)(23)(34)(45)(51). 

In summary, we have explained how the evaluation of the correlator G n+ 2 is reduced to 
algebraic manipulations by the insertion procedure combined with superconformal symmetry 
and harmonic analyticity. As in the four-point case, no integral needs to be done once the T 
block is known. Nevertheless, the amount of algebra is fairly large, so that we have resorted to a 
Mathematica script. It remains to restore the integrations over the isnertion points xq,xq' with 
the IR-modified measure. The scalar (parity-even) part of the result for G 5 has the concise form 
displayed in Eq. (6.7). 

A special comment is due here on the pseudo-scalar (parity-odd) part of the correlator. Our 
calculation of G5+2 does indeed produce such a (rather complicated) part. But this does not mean 
that G5 will have a parity-odd part. It must drop out after the integration over the insertion 
points Xq, xq>. The explanation is given in Appendix A. 5. 

A. 5 Parity properties of the scalars in the J\f = 2 theory 

Here, following Ref. [47], we argue that there exists a parity assignment for the fields of the 
M = 2 vector and hypermultiplets, such that the hypermultiplet scalars are true scalars (not 
pseudo). With this assignment, all our operators made of hypermultiplets are scalars, and their 
correlators should not contain a parity-odd part. 

The components of the two multiplets are contained in the G-analytic superfields q + , Eq. (A. 8), 
and V ++ , Eq. (A. 13). The parity assignments of Ref. [47] for the superspace coordinates are 

P: x' =x , x' = -x, (0+)' = 0+*, (6+)' = -e + \ (ut)' = ut, (AM) 

while the superfields remain inert, 

q + '(x',e',u') = q + (x,e,u), V ++ '(x',$',u') = V ++ (x,e,u). (A.62) 

In terms of the bosonic physical fields, these assignments imply that the hypermultiplet scalars 
(p 1 are true scalars, the gluon A^ is a polar vector, while the complex vector multiplet scalar ip is 
a mixture of a true and a pseudo-scalar. In addition, the fermion fields transform in an unusual 
way: the hyperinos transform into each other (up to signs), ip -H- K, and the gluinos A* -H- A*. The 
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latter relation means that, e.g., A 1 -H- A 1 = A2 = (A 2 )*, which again differs from the traditional 
assignment A* -B- (A*)*. 

These rules can be tested for consistency by inspecting the Yukawa couplings in the M = 4 
Lagrangian. The gauge/matter coupling reads 

f dud 4 xd 4 9 + Ti(V ++ [q + ,q + }) =► Tr[0 4 ({/t, A;} + {$, A,})] + c.c. (A.63) 

We see that the above assignments allow 0* to stay inert under parity. At the same time, the 
Yukawa coupling from the gauge sector is 

[ d 4 xd 4 9 Tr(W 2 ) =► f d 4 xTr(p{\ i ,\ i } + cp{\ t ,\ t }) . (A.64) 

Here the combination of fermions accompanying each boson is chiral, therefore </?■<->■ (p. 

The above parity assignments mean that the hypermultiplet composite operators of the type 
O = Tr(q + ) k \ e=0 = Tr(0 ll (x) . . . 0* fc (x))w^ . . .uf are all true scalars. Thus, the correlators 
(OOOO) that we are considering cannot have a parity-odd part. 

The reason why we see such a part in the correlators G n+ \ and G n+ 2 was explained after 
Eq. (4.9). It is due to the insertion of the complex (chiral) form of the SYM Lagrangian. In 
it we find pseudo-scalar terms, for example iFF, which are responsible for the parity-odd part 
in the correlators with insertions. But at the final stage of the calculation, the integration over 
the insertion point will eliminate all such terms, which are total derivatives. Indeed, we have 
already encountered this phenomenon in Sect. 4. The correlator with one insertion G n+ \ has the 
form (4.9), where we clearly see a pseudo-scalar in the last line. However, it drops out after the 
integration over the insertion point. 

B Four-point correlators of operators of weight k 

Let us consider the four-point correlator of protected half-BPS operators of weight k. At the 
lowest level of the 9 expansion they are built from k elementary scalar fields, Tr(0 fc ). Such 
correlators are discussed in detail in [44, 45], using the method of [23]. We consider the following 
M = 2 hypermultiplet projection 

G A]k = (O k {l)d k {2)O k (3)O k {4)) (B.l) 

where in this appendix O k = 2 fc//2 Tr(g + ) fc . 

The (connected, planar) tree-level correlator is given by (see Fig. 9) 



fc-i 



G<g = k 4 N 2 c k ~ 2 Y, X m Y k -' m , (B.2) 



m=l 



where the notation was introduced for the harmonic and space-time propagator factors 

x _ (12) (34) ^ y _ (14) (32) ^ (R3) 

{2ir) x 12 x 34 (27TJ x 14 x 2 z 
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Figure 9: Tree graphs for operators of weight k. 
The loop corrections have the general factorized form 

fe-2 

G (ioo P s) = n , j2 f k m ( x )X m Y k ~ m ~ 2 

with the polynomial prefactor 

-R! = sX 2 + (t - s - 1)XY + Y 2 
involving the conformal cross-ratios 



2 2 ' 
x 14 x 23 



x 13 x 24 
x 14 x 23 



The coefficient functions J-(x) will be specified below. 

The result of [44] for the one-loop correction in the planar limit is 



fc-i 



G^ = k 4 N 2k - 2 n' 2ax 2 14 x 2 23 g(1234)J2 xm ~ lYk ~ m ~ 1 



m=l 



(B-4) 



(B.5) 



(B.6) 



(B.7) 



where the one-loop box integral <?(1234) is defined in (5.9). We see that in this case the general 
amplitude (B.4) becomes completely degenerate, with all J-^Xx) oc Xi 4 x| 3 y(1234). Going to the 
light-cone limit x 



J/oq — X 



12 — ^23 ~~ ^34 — ^41 



Xai — > 0, we find that s remains finite while t — > oo, hence 



K' -»• tXY . 

As a consequence, Eq. (B.7) simplifies to 

fe-i 
G-W = fc 4 jV c 2fc " 2 2az? 3 a& #(1234) ^ X m r fc - m + subleading terms , 



(B.8) 



(B.9) 



m=l 



Dividing this expression by the tree-level correlator (B.2), we obtain the same ratio as in the case 
k — 2 at one loop, see (5.18). Thus, the ratio does not depend on the value of k, up to one loop. 
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The same pattern is found at two loops. According to [45] , the two-loop coefficient functions 
in (B.4) are 

9 4 

^ m = /o_-2\2 I [ ( -'m x i3 x 24 + (2sC m — C m )x u x 23 + (2C m j s — C m ) x 12 x 34 \ [#(1234)] 

+ (Cm - Ci)2[x 2 13 h(123; 134) + 4^(124; 234)] 
+ (Ci - C- 1 )2[x? 4 /i(124; 134) + 4^(123; 234)] 
+ (C£ - C^)2[x? 2 /i(123; 124) + 4^(134; 234)]} , (B.10) 

with the two-loop integral h defined in (5.9) and with color factors C given in [45]. Going to the 
light-cone limit, only the first term from the first line and the term from the second line survive. 
Further, in the planar limit the remaining color factors simplify to (see [45]) 

C d m = k*N 2k , C c m = 2k 4 N 2k . (B.ll) 

As a result, we find that on the light cone all the relevant two-loop coefficients become equal, 

J*{x) = a 2 k*N 2k ~ 2 {(x? 3 ^ 4 ) 2 [^(1234)] 2 + 2x 2 13 x 2 24 [x 2 13 h(123; 134) + o: 2 4 /i(124; 234)]} . (B.12) 

Substituting this result into (B.4) and taking account of (B.8), we obtain the two-loop correction 

Cg ={(x 2 l3 x 2 24 ) 2 [g(123A)] 2 + 2x 2 l3 x 2 2A [x 2 13 h(123; 134) + x 2 2A h(124; 234)]} 
fe-i 
k i N 2k-2 J2 x m Y k - m + subleading terms . (B.13) 



x a 

m=l 



Finally, dividing by the tree-level correlator (B.2), we obtain the same ratio as in the case k = 2 
at two loops, see (5.18). This confirms that the ratio does not depend on the value of k up to 
two loops. 

In conclusion, we can claim that the duality relation 

lim G A]k /G { ;i=(A A]k /A^\0^) (B.14) 

holds for any weight k of the half-BPS operators. 

B.0.1 Relation to Wilson loops 

From our analysis it follows that the correlator in the planar limit has the following universal 
form on the light-cone 

(g fel (l)^ 2 (2)g fe3 (3)g fc4 (4)) = (g fel (l)^ 2 (2)g fe3 (3)g fc4 (4))(°) [W(x 1 ,x 2 ,x 3 ,x 4 )} 2 , (B.15) 

where W(xi) is a light-like Wilson loop in the fundamental representation of SU(N C ) evaluated 
along the contour [xi,^] U [£2,2:3] U [^3,^4] U [x4,Xi], 

Wfa) = — (0|TrPexp (ig j> dx ■ A(x) J |0) . (B.16) 
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iY, 



Efc; 



Figure 10: In the planar limit each propagator is replaced by a pair of arrowed lines and the vertex is 
replaced as shown in Fig. 11. Each line with an arrow corresponds to a Wilson line. 





Figure 11: Each vertex q in the planar limit is replaced by k pairs of arrowed lines. 

The diagrammatic derivation of the above relation is shown in Figure 10. As explained in 
the parallel paper [28], each propagator connecting a pair of adjacent points Xi and Xi + ± is 
approximated by a free propagator multiplied by a Wilson line along the segment [xi,Xj + i], 
evaluated in the adjoint representation, VFadj^^i+i]- Then, the vertex at point Xi contains ki 
Wilson lines with their color indices contracted to ensure that the total color charge is zero. The 
Wilson lines in the adjoint and in the fundamental representations are related to each other as 



(W 3di [x i ,Xi +1 ]) ab t b = W iund [xi, x i+1 ]t a W iund [x i+1 , Xi] 



(B.17) 



or equivalently (for the gauge group U(N)) 

(Wf U nd[^,^+i])ij(VFf und [xi + i,Xi]) M = (t a ) kj (W ad j[xi,Xi +1 ]) ab (t b )ii . (B.18) 

In the multi-color limit, we can use the last identity to replace a Wilson line in the adjoint 
representation by a pair of two parallel fundamental Wilson lines with opposite orientations. 
This is shown in the middle panel of Fig. 10. In this way, we obtain a collection of closed 
cycles. We observe that all cycles but two have a backtrack shape, i.e. the corresponding contour 
encircles a zero area. We denote such contour by C U C~ x . Only two cycles go through all cusp 
points Xi with different orientations. Notice that the Wilson lines satisfy the unitarity condition 



Wc(W ( 



c 



it 



WcWc-i = 1 • 



(B.19) 



As a consequence, each backtrack cycle reduces to 1 and we arrive at the right-hand side panel 
in Fig. 10. It contains only two cycles, each corresponding to a Wilson loop in the fundamental 
representation. 
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